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Chapter  1 

INTRODUCTION  AND  SUMMARY 


1-1.  Introduction 

This  Volume  III  of  this  final  technical  report  of  four  volumes 
documents  a  surface  integration  radome  analysis  computer  program  written 
in  Fortran  IV  for  use  on  the  Cyber  70/74  computing  system  at  Georgia 
Institute  of  Technology  and  the  IBM  3033  computing  system  at  Johns  Hopkins 
University  Applied  Physics  Laboratory.  The  program  was  developed  at 
Georgia  Institute  of  Technology  over  the  past  three  years  under  grant 
AFOSR-77-3469  and  documented  herein  under  the  cognizance  of  R.  C.  Mallaleiu 
(APL  Contract  60153) . 

The  analysis  package  described  was  used  during  the  research  to 
analyze  the  antennas  and  radomes  as  described  in  Volumes  I  and  IV.  Its 
documentation  was  done  in  conjunction  with  an  on-going  radome  technology 
program  at  JHU/APL.  It  is  intended  to  serve  as  part  of  a  technology  base 
for  the  radome  technical  community. 

This  report  is  organized  by  chapters,  where  each  chapter  describes 
the  main  program  or  one  subprogram  not  already  described  in  Reference  1 
(Ray  Tracing  Formulation) .  The  main  program  (Chapter  2)  described  herein 
differs  only  slightly  from  that  in  Reference  1.  Only  six  new  subroutines 
are  required  for  the  surface  integration  formulation  as  described  in 
Chapters  3-8.  References  cited  in  each  chapter  are  listed  therein.  Each 
chapter  is  terminated  with  the  program  listing. 

This  software  us  currently  being  used  in  a  parametric  investigation 
of  radome  analysis  methods,  and  additional  information  concerning  its 
speed  and  accuracy  is  presented  in  Volume  I  [2). 
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1-2.  Description  of  the  Analysis 

The  basis  of  analysis  is  illustrated  in  Figure  1-1.  The  inner  sur¬ 
face  of  the  radome  is  chosen  as  one  surface  of  integration  in  the 
Lorentz  reciprocity  integral  (upper  left  in  figure) .  A  surface  enclos¬ 
ing  the  antenna  and  extending  into  its  interior,  as  illustrated,  comprises 
the  second  surface.  Together,  +  S2  enclose  the  source-free  Volume  V 
indicated;  hence,  the  Lorentz  surface  integral  is  identically  zero  and  the 
integral  over  equals  the  negative  of  the  integral  over  S^. 

Consider  the  surface  S2  more  closely.  The  surface  integral  over 
is  zero  except  over  that  part  of  S2  which  is  placed  across  the  waveguiding 
structure  that  connects  the  Source  "a"  to  the  radiating  (flared)  part  of 
the  antenna.  Call  this  surface  S^.  If  there  can  be  defined  a  single 
dominant  mode  in  the  waveguide  when  Source  "a"  is  activated,  then  voltage 
and  current  V  ,  I  can  be  defined  at  this  terminal  plane  [4] .  When  Source 

3  3. 

"b"  is  activated,  voltage  and  current  V  ,  I  at  S'  can  also  be  defined; 

b  b  2 

in  fact,  V,  is  the  "received  voltage”.  The  received  current  I,  is  related 
b  b 

to  V  by  a  linear  impedance  relationship 


(1) 


where  Z  is  the  impedance  seen  at  S'  looking  toward  Source  "a"  (sinusoidal 
steady  state  assumed;  time  variations  of  the  form  understood  and 

suppressed) .  Also,  I  and  V  are  related  by 


V 

a 


(2) 


2 


THEORETICAL  BASIS  OF  RADOME  ANALYSIS. 


where  Z^  is  the  impedance  seen  at  looking  to  the  right  in  Figure  1-1. 

Combining  these  results  yields  the  desired  expression  for  the  received 

voltage  V  ;  viz., 
b 


V  =  V 
b  REC 


Z.Z 
1  a 


V  (Z  +  z  ) 
a  1  a 


(E  X  H 

—a  — b 


E,  x  H  )  •  n  da 

— b  —a 


(3) 


Note  that  the  unit  normal  is  directed  positively  outward  from  volume  V 
as  dictated  by  Gauss’  theorem. 

When  Source  "b"  in  Figure  1-1  is  removed  a  great  distance  from  the 
antenna/radome  structure,  the  fields  of  "b"  approach  those  of  an  electro¬ 
magnetic  plane  wave  (target  return) .  The  practical  analysis  approach  then 
takes  the  form  shown  in  Figure  1-2.  The  inner  radome  surface  is  divided 
into  a  number  of  contiguous  elemental  areas  AA,  each  of  which  is  repre¬ 
sented  by  a  sample  point  P'  at  its  center.  The  fields  at  P'  are 

assumed  to  be  those  present  there  in  the  absence  of  the  radome  and  are 
found  by  aperture  .integration,  the  theoretical  basis  of  which  is  the  Huygens- 
Fresnel  principle  [5]  as  stated  by  the  lower  integrals  in  Figure  1-1.  The 

fields  E  ,  H  at  P'  are  found  by  applying  the  normal  voltage  transmission 
R  R 

coefficients  [6]  to  the  plane  wave  incident  on  the  outside  at  point  P. 

The  received  voltage  is  found  by  summing  all  the  contributions  as  indicated 
in  Equation  (3 ) . 

The  method  of  analysis  indicated  by  Equation  (3)  is  exact;  however, 

certain  approximations  are  necessarily  introduced  in  its  implementation. 

The  fields  E H  should  correctly  include  reflections  from  the  inner 
-T  -T 

radome  surface.  The  use  of  the  flat  panel  transmission  coefficients  to 
transform  the  incident  plane  wave  at  P  to  P'  is  an  approximate  method 
based  on  the  theory  of  geometrical  optics  (zero  wavelength)  and  whose 
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ac'  ^ucy  deponds  upon  the  radius  of 
accuracy  of  the  method  also  depends 
represent  the  radiating  aperture  as 
computational  speed  of  the  analysis 
of  these  samples. 
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Chapter  2 


PROGRAM  SIIRACP 

2-1.  Purpose:  SIIRACP  is  a  Fortran  computer  program  used  to  analyze 

the  effects  of  a  tangent  ogive  radome  on  the  performance  of  a  mono- 
pulse  aperture  antenna.  It  consists  of  a  main  program  and  28  sub¬ 
routines,  22  of  which  are  identical  to  those  used  in  Program  RTFRACP 
[1] .  It  uses  complex  arithmetic  and  requires  66600  octal  words  of 
core  memory  for  execution  on  the  CDC  Cyber  70  system  (60-bit  words) 
at  Georgia  Institute  of  Technology.  Execution  time  to  compute  bore- 
sight  error  on  the  Cyber  70  is  approximately  255  seconds  per  look 
direction  when  the  small  antenna  aperture  is  represented  by  7  x  7  = 

49  sample  data  points  and  the  radome  is  represented  by  826  sample 
points;  i.e.,  approximately  1.26  millisecond  per  aperture  sample 
point  per  radome  sample  point. 

The  computer-aided  radome  analysis  uses  a  receiving  formulation 
based  on  the  Lorentz  reciprocity  theorem  as  described  earlier  [1,2].  The 
voltage  produced  at  the  terminals  of  a  linear  antenna  by  an  incident  plane 
wave  is  given  by 

„  re 

VR(k)  =  jj  <£r  x  £r  '  -r  x  V  •  n  da  (1) 

s 

where  E^,  are  the  fields  produced  on  the  surface  S  enclosing  the  antenna 
when  the  antenna  is  transmitting;  E  ,  H  are  the  incident  fields  produced 
on  S  by  the  incident  plane  wave  or  perturbations  thereof;  k  is  a  unit  vec¬ 
tor  which  points  from  the  antenna  toward  the  direction  from  which  the  plane 


wave  arrives;  and  n  is  a  unit  vector  normal  to  the  surface  S  and  pointing 

★ 

into  the  source-free  region.  The  fields  E^,  are  taken  to  be  those  pro¬ 
duced  by  the  planar  aperture  on  surface  S  when  the  antenna  is  transmitting 
in  the  absence  of  the  radome.  The  geometrical  optics  approximation 


(2) 


is  used  to  generate  the  magnetic  field  in  the  aperture  from  the  aperture 
illumination  specified  by  E^. 

The  surface  S  is  taken  to  be  the  inner  surface  of  the  radome.  At 

each  sample  point  P'on  this  surface,  the  plane  wave  fields  E  ,  H  incident 

R  R 

from  the  outside  are  weighted  with  the  flat  panel  normal  voltage  transmis¬ 
sion  coefficients  as  determined  by  the  radome  wall  configuration,  the  angle 
of  incidence,  and  the  plane  of  incidence.  The  fields  E^,  a't  P'  are  found 
by  aperture  integration.  The  individual  contributions  are  summed  up  as 
indicated  in  Equation  (1)  and  was  illustrated  in  Figure  1-2. 

The  parameters  of  the  tangent  ogive  radome  are  indicated  in  Figure 

2-1.  The  outside  base  diameter  D  and  fineness  ratio  F  determine  the 

os  os 

outside  length  according  to 


F  =  L  /D 
os  os  os 


A  similar  relation  holds  for  the  inside  dimensions;  viz., 


F  =  L.  /» 
is  is  is 


(3) 


(4) 


By  choosing  h  this  way,  the  minus  sign  in  Figure  (1-1)  is  removed. 


A  .  . 


The  radius  of  curvature  of  the  outside  wall  R  is  given  by 

os 

R  =  F  D  /  sin  ( tt  -  2  Tan  1(2F  ))  (5) 

os  os  os  os 

and  the  dimension  B  is  given  by 

B  =  R  -  D  /2  (6) 

03  OS 

The  placements  of  a  bulkhead  (bottom  disk)  and  metal  tip  (top  disk)  can 
be  specified  by  and  Z  ,  respectively.  The  thickness,  dielectric 

constant,  and  loss  tangent  of  the  wall  may  also  be  specified  for  up  to 
N=5  layers.  The  radome  is  assumed  to  be  a  body  of  revolution  with  uni¬ 
form  wall  dimensions  independent  of  location.  The  dashed  cylindrical 
shape  of  a  diameter  in  Figure  2-1  was  used  earlier  to  simulate  a  laser- 
induced  defect  and  is  not  pertinent  here. 

The  subroutine  which  generates  the  antenna  aperture  fields  represents 
throe  typos  of  antennas:  circular  or  square  aperture  with  tapered  (cos  x) 
illumination  and  any  one  of  four  polarizations  (vertical,  horizontal,  RHC, 
LHC);  flat  plate  antenna  with  tapered  illumination  and  vertical  polariza¬ 
tion.  For  either  antenna,  the  fields  are  computed  for  one  of  three  selected 
channels:  sum,  azimuth  difference,  elevation  difference.  Inputs  include 

the  number  of  samples  N  ,  N  and  the  aperture  diameter  D  /X  in  wavelengths. 

X  Y  AP 

The  antenna/radome  orientation  is  specified  according  to  the  para¬ 
meters  defined  in  Figure  2-2.  The  angle  $  selects  the  plane  of  scan  of 
the  radome  tip  with  respect  to  the  antenna  coordinate  system:  <p  =  0° 
selects  the  azimuth  plane;  =  90°  selects  the  elevation  plane.  The 
angle  0  scans  the  tip  in  the  selected  plane. 

ii 
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The  program  computes  boresight  errors  in  the  azimuth  and  elevation 

1  lanes  of  the'  antenna.  The  radome  orientation  is  specified  by  f>  and  if  . 

P  L 

The  first  tarqet  return  (plum-  wave)  is  made  to  arrive  from  the  direction 

k  -  x  sin  *i  t  y  sin  i)  +  /.  tl  -  2  :;in‘  f)  (7) 

1  A  os  A  os  A 


where  •'  is  tin 

os 

l  n  i  t  i  a  ) 

Specified  offset 

angle ; 

<-•  •  g  • 

,  2°.  The  voltage 

received  by  each 

channel 

is  computed  and  : 

stored . 

The 

second  return  is 

made  to  arrive  fr 

i  un 

k  x  I -sin  ■  )  +  v  (-sin  "  )  t  7,  r’l  -  sm‘  6  (8) 

..’A  s  A  i  is  A  os 


and  the  voltages  art'  again  computed.  Ttie  data  from  those  two  points  are- 
used  to  construct  a  linear  tracking  model  in  the  two  [lanes,  and  a  direc¬ 
tion  of  ariival  k  l  .  |  r»*di  -ted  which  will  yield  null  indications  in  both 


;  1  dfi'- 

.  Tin  i-rric, 

•ss  is  r.  i  eat  ed 

tint  i  1  a  i It  s 

i  red 
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r  a  maximum 
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t ions  i s  .  x 

eeedl 

■d .  I'pon  completion 

,  tin* 

<  ut :  »it 

k  i ndicat •  : 

:  t.  lie  dll  eC  t  1  >n 

trim  wli  i  eh 

t  h.  [ 

Ian*  arrives  which 

yie] ds 

i  r  i  •  1  • 

>  t  i  :  i  .i  1  lit' 

•;  i  •:!.*  i  nd  1  ca  t  I  on 

I t  i  Hid 

!•  r. 

■present  the>  bore-sight 
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Options  are  also  provided  whereby  principal  plane  patterns  as  shown 
in  Figure  2-3  and  additional  outputs  around  boresight  can  be  computed  and 
printed.  These  options  are  useful  when  preparing  software  for  a  new  type 
of  antenna  and  to  ensure  correct  operation  whenever  curious  results  are 


obtained. 

2-2.  Usage:  Line  No. 

DATA  APIN/O ./  47 

DATA  ZBOTIN/O . 00/  49 

DATA  RADIUS/ 1E0/  52 

DATA  THETAA,  PHIA,  AGAM3A/0 .0, 90 . 0, 0 . 0/  53 

DATA  NX,  NY,  NYE,  NXY/4 ,4 , 1 , 512/  56 

DATA  MY/1/ , NREC/61  57 

READ  (5,6)  TITLE  63 

READ  (5,*)  GRAF3D,  GRAFSA,  GRAFTR,  GRAFRV,  SUPPRS,  IPENCD, 

SQUARE  65 

READ  (5,*)  NFINE,  NPHI ,  NTHE,  DIAOS,  RAIN,  RRIN ,  ZTOPIN,  FREQ, 
OSANG  67 

READ  (5,*)  LMAX,  DMRAD,  IOPT,  RAPMAX,  VAIRM,  IPOL,  ICASE, 

N,  IPWR,  UMAX,  NXE  76 

READ  (5,*)  DSTHIN,  DSPHIN,  NTHMIN,  NPHI MIN  79 

READ  (5,*)  DIN ( I) ,  ER(I),  TD(I)  (1=1, N)  115 

READ  (5,*)  FINR(I)  (1=1, NFINE)  124 

READ  (5,*)  PHI (I )  (1=1,  NPHI)  127 

READ  (5,*)  TOETA(I)  (1=1, NTHE)  129 
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£  =  elevation  component 

A 

a  =  azimuth  component 


2-3.  Arguments 


a.  Inputs.  Units  of  arguments  on  input  are  distances  in  inches, 
angles  in  degrees,  and  frequency  in  gigahertz,  unless  otherwise  noted. 

Units  of  arguments  passed  to  subroutines  are  centimeters,  radians,  and 
gigahertz.  An  asterisk  is  used  to  denote  those  DATA  arguments  that  do 
not  normally  need  to  be  changed  by  the  user. 

APIN*  -  Height  of  a  cylindrical  base  section  of  the  tan¬ 

gent  ogive  radome .  It  is  no  longer  included  in 
the  ray  tracing  algorithms  and  should  not  be 
changed  from  its  zero  value. 

ZBOTIN  -  Distance  from  base  of  tangent  ogive  radome  to 

missile  bulkhead  (Figure  2-1) . 

RADIUS*  -  The  radius  R  used  in  the  far  field  factor  e  ^kR/R 

by  Subroutine  FAR.  Do  not  change. 

THETAA*  -  Angle  0  between  z-axis  and  the  position  vector  r 

a  -a 

to  the  antenna  origin.  This  angle  was  used  in 

earlier  work  to  locate  the  antenna  origin  in  the 

reference  system  using  spherical  coordinates 

(r  ,  8  ,  $  ) .  Do  not  change, 
a  a  a 

PHIA*  -  Angle  <f>  between  the  projection  of  z  axis  onto 

3  A 

the  xy-plane  and  the  x-axis.  Do  not  change. 

AGAM3A*  -  Angle  between  z^-axis  and  z-axis  in  Figure  2-2. 

Do  not  change. 

NX, NY  -  Integer  powers  of  two  equal  to  the  number  of  sample 

points  in  the  antenna  aperture;  e.g.,  16,  32,  64, 
etc.  Changing  NX  and  NY  necessitates  compatible 
changes  in  Lines  16-18. 


IS 


NYE 

(NXE) 


Integer  powers  of  two  which  specify  the  expanded 


NXY 


NREC 


NS 


MX, MY 


number  of  sample  points  desired  when  computing 
the  transmitting  patterns  of  the  antenna  by 
inverse  Fourier  transforming  the  aperture  fields. 
Subroutine  JOYFFT  provides  this  capability  of  in¬ 
creased  resolution  in  one  or  both  dimensions. 

Changes  in  NXE,  NYE  necessitate  compatible  changes 
in  Lines  16,  20,  22  and  23.  Note  that  NXE*NYE<NX*NY 
and  either  NXE<NX  or  NYE<NY . 

Integer  power  of  two  used  by  Subroutine  JOYFFT  for 
dimension  of  complex  working  array  XYFFT .  Note 
that  MX*NX<NXY  and  MY*NY<NSY .  See  below  for  MX 
and  MY . 

Integer  variable  equal  to  the  number  of  points 

at  which  to  compute  the  receiving  pattern  in  either 

principal  plane.  The  received  voltage  is  computed 

at  points  6 ^  equally  spaced  in  sinP,  where  9  is 

the  angle  measured  from  the  z^-axis  as  indicated 

in  Figure  2-3,  vhere  sin  0  =  -KMAX  +  ( 1-1) *2*KMAX/ 

NREC,  and  where  KMAX  =  sin  0  <1.0. 

max 

Not  used.  It  was  originally  used  by  Subroutine 
RECBS  .  Do  not  remove . 

Integer  powers  of  two  equal  to  the  magnification 

factors  desired  in  the  k  (H-i lane)  and  k  (E-plane) 

x  y 

directions,  respective ly ,  of  tin-  transmitting  an¬ 
tenna  patterns.  Note  that  the  restrictions  MX*NY<NXY 
and  MY*NY<NXY  must  be  observed.  The  data  cited 


16 


above  indicates  increased  resolutions  in  the  NX 


TITLE 


GRAF3D 


GRAFSA 


GRAFT R 


GRAFRV 


direction  of  MX=16  and  no  magnification  (MY=1)  in 
the  NY  direction.  Consequently,  note  that 
NXE=MX*NX=256 . 

A  Hollerith  string  of  up  to  72  characters  which 
describes  briefly  the  analysis  being  done.  A 
format  of  18A4  is  specified  and  should  work  for 
machines  with  word  length  greater  than  or  equal 
to  32  bits.  The  dimension  of  TITLE  (Line  31) 
should  be  at  least  18 . 

-  A  logical  variable  used  to  control  the  plotting 

of  the  incident  fields  on  the  antenna  aperture. 

This  feature  has  been  removed  from  the  program, 

and  GRAF3D  should  always  be  FALSE. 

A  logical  variable  which  (if  TRUE)  controls  the 

plotting  of  the  transmitting  power  patterns  of 

the  antenna  as  follows:  E-plane  sum.  E-plane 

difference  equation  (A  )  ,  H-plane  sum,  and 

EL 

H-plane  difference  azimuth  (A  ) .  The  radome 

AZ 

is  absent. 

A  logical  variable  which  controls  the  plotting  of 
the  amplitude  and  phase  of  the  antenna  aperture 
fields  in  the  following  order: 

Exi'  ey£'  exael'  eyael'  exaaz'  eyaaz” 

A  logical  variable  which  controls  the  plotting  of 
the  receiving  patterns  of  the  antenna  with  radome 
in  the  same  order  as  specified  under  GRAFSA  above 
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SUPPRS 


A  logical  variable  which  controls  the  printing  of 


IPENCD 


SQUARE 


NFINE 


NPHI 


NTHE 


DIAOS 


numerous  results.  When  TRUE,  the  printing  of  these 
numerous  results  are  suppressed.  This  feature  is 
convenient  to  aid  in  debugging  new  portions  of 
software  prior  to  making  production  runs. 

An  integer  variable  which  selects  pen  and  paper 
for  the  Cal  comp.  This  variable  may  be  system 
dependent.  For  the  Cyber  70,  IPENCD=00  yields 
ballpoint  pen  and  11”  wide  plain  paper;  IPENCD=40 
yields  a  heavier  ink  pen  and  the  same  paper. 

Logical  input  variable  which  selects  a  square 
aperture  (TRUE)  in  Subroutine  TRECNF. 

Integer  variable  equal  to  the  number  of  fineness 
ratios  to  be  considered  for  the  tangent  ogive 
radome;  e.g.,  NFINE=1. 

Integer  variable  equal  to  the  number  of  scan  planes ; 
e.g.,  NPHI=2 . 

Integer  variable  equal  to  the  number  of  angles  in 
each  scan  plane  at  which  to  compute  boresight 
errors,  etc.  Note:  The  program  is  set  up  to 
iterate  on  fineness  ratio,  scan  plane,  and  scan 
angle  as  outer  loop,  middle  loop,  and  inner  loop, 
respectively.  Therefore,  for  each  of  NFINE  fine¬ 
ness  ratios,  the  analysis  will  be  done  for  NTHE 
scan  angles  in  NPHI  different  scan  planes. 

Real  variable  equal  to  the  outside  base  diameter 
(in.)  of  the  radome.  See  Figure  2-1. 
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RAIN 


Real  variable  equal  to  the  distance  (in.)  from 


RRIN 


ZTOPIN 


FREQ 


OSANG 


LMAX 


DMRAD 


I  OPT 


RAPMAX 


VAIRM 


the  gimbal  point  to  the  antenna  aperture. 

Real  variable  equal  to  the  distance  (in.)  from  the 
gimbal  point  to  the  base  of  the  radome . 

Real  variable  equal  to  the  distance  (in.)  from 
the  base  of  the  radome  to  the  face  of  a  metal  tip 
on  the  radome . 

Real  variable  equal  to  the  frequency  of  operation 
in  gigahertz. 

keal  variable  equal  to  the  offset  angle  in  degrees 
at  which  the  first  target  return  is  to  arrive  on 
the  antenna;  e.g.,  OSANG=3.0. 

Integer  variable  equal  to  the  maximum  number  of 
iterations  allowed  by  Subroutine  RECBS  in  com¬ 
puting  boresight  error;  e.g.,  LMAX=5. 

Real  variable  equal  to  the  tolerance  in  milliradians 
allowed  on  computing  boresight  error;  e.g.,  DMRAD=0 . 1 . 
Integer  variable  which  selects  the  polarization 
of  the  incident  plane  wave  as  follows: 

1.  Linear,  elevation  component 

2.  Linear,  azimuth  component 

3.  Right  hand  circular 

4.  Left  hand  circular 

Real  variable  equal  to  the  maximum  radius  (in.) 
of  the  antenna  aperture. 

Real  variable  equal  to  the  maximum  amplitude  of 
sum  channel  received  voltage  without  radome.  Any 
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real  value  can  be  entered  for  this  variable  since 


I  POL 

I  CASE 


N 


IPWR 


KMAX 

NXE 

DSTHIN, 
DS  PH  I N 


a  subsequent  program  modification  (Lines  J43-3L2) 
causes  VAIRM  to  be  computed  automatically. 

Integer  variable  which  selects  the  polarization  of 
the  antenna  when  ICASE=1  according  to  the  same  code 
as  used  above  for  IOPT. 

Integer  variable  which  selects  the  type  of  antenna 
aperture  for  the  analysis:  I CASE- 1  or  selects 
a  circular  or  square  aperture  with  tapered  illumi¬ 
nation:  ICASE=3  selects  a  flat  plate  antenna  with 
programmed  illumination.  See  Subroutine-  TRECNK  in 
Chapt:  r  4. 

Integer  variable-  equal  to  the  number  of  layers 
(up  to  5)  in  the  radome  wall.  For  cases  where 
more  than  h  layers  are  required,  the  dimensional 
arrays  on  Line  37  must  be  changed  to  NN=N+1. 

Integer  variable  which  selects  the  component  for 
which  to  compute  the-  transmitting  power  i  attorns 
as  follows: 

1 .  Elevation  Comjxments 

2.  Azimuth  Component 

3.  Total  power 

Real  variable  equal  to  the  sine  of  the  maximum 
angle  at  which  receiving  patterns  are  to  he  computed. 
Integer  variable  used  by  JOYFFT  as  explained  above. 

Rt al  variables  equal  to  the  distance  between  adja¬ 
cent  sample  points  on  the  radome  surface  in  the 
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longitudinal  (6)  and  circumferential  (<*.)  di¬ 
rections,  respectively.  See  Chapter  3. 

NTHMIN,  -  Integer  variables  equal  to  the  minimum  acceptable 

NPHIMIN 

number  of  radome  sample  points  in  the  two  direc¬ 
tions  . 

DIN,ER,TD  -  Subscripted  real  variables  equal  to  the  thickness 
(in.),  dielectric  constant  (e^),  and  loss  tangent 
(tan  <$)  of  each  layer  of  the  radome  wall.  1=1 
corresponds  to  the  first  layer  and  is  the  layer  on 
exit  side  of  the  wall.  Layer  N  is  the  first  layer 
encountered  by  the  incident  plane  wave.  See  Sub¬ 
routine  WALL. 

FINR  -  Subscripted  real  variable  equal  to  NPINE  fine¬ 

ness  ratios. 

PHI  -  Subscripted  real  variable  equal  to  NPHI  angles 

(degrees)  which  specify  the  scan  planes. 

THETA  -  Subscripted  real  variable  equal  to  NTHE  angles 

(degrees)  which  specify  the  scan  angles  in  the 
scan  plane. 

b.  Outputs.  The  parameters  of  analysis  which  are  computed  and 
outputted  by  the  program  depend  on  whether  SUPPRS  is  true.  In  what 
follows,  it  is  assumed  that  SUPPRS=FALSE  so  that  all  possible  outputs  are 
obtained.  Since  many  of  the  original  input  parameters  arc  printed  di¬ 
rectly,  only  those  parameters  not  already  explained  above  will  be  included 
below.  Additional  clarification  may  be  found  in  Section  2-6. 

TABLE  -  Logical  variable  which,  if  TRUE,  causes  a  look-up 

table  to  be  used  in  computing  transmission  coeffi¬ 
cients.  When  SUPPRS=FALSE,  an  abbreviated  table 
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ANGLE 


TPERI , TPARI  - 


RPERI , RPARI 


KXMAX 


DXWL 
KXM , KYM 


of  transmission  coeffi cients  of  the  radome  wall 
is  printed  by  Subroutine  WALL  with  variables  as 
explained  immediately  below. 

Real  variable  equal  to  the  angle  of  incidence 
(degrees)  of  the  plane  wave  on  a  plane  sheet  of 
infinite  extent  having  the  layered  configuration 
specified  for  the  radome  wall.  The  entries  in 
the  table  are  computed  at  250  equal  increments 


in  sin  0.,  but  only  every  fifth  result  is  printed. 


Complex  variables  equal  to  the  normal  voltage  trans¬ 
mission  coefficients  of  the  sheet  for  the  two  cases 


of  EL  perpendicular  to  the  plane  of  incidence  (T^) 


and  E.  parallel  to  the  plane  of  incidence  (T..) 


— 1 


In  the  printed  table,  the  power  transmission  coef¬ 


ficients  !t±  I  are  |T,,|  *"  ate  printed;  adjacent  to 


each,  the  phases  of  and  T|j  are  also  printed. 


Complex  variables  equul  to  the  reflection  coeffi¬ 


cients  R^  ,  R|.  of  the  plane  dielectric  sheet. 


2  2 

Actually,  |R  and  ! R ..  I  are  printed,  accompanied 


by  the  nhases  of  R  and  R... 

•  l  i 


Real  variable  equal  to  the  folding  wavenumber  asso¬ 
ciated  with  sampling  the  aperture  fields  according 
to  KXMAX  =  l./(2(Ax/A),  where  x  is  the  distance 
between  samples. 

Real  variable  equal  to  Ax/A. 

Real  variables  equal  to  the  folding  wavenumbers  of 
the  principal  plane  patterns  after  magnification 


1 


2  2 


for  increased  resolution.  KXM=KYMAX*NXE/ ( MX*NX) 
and  applies  to  the  H-plane. 

KYM=KYMAX*NYE/(MY*NY)  and  applies  to  the  E-plane. 
Usually,  the  expanded  dimension  NXE  and  magnigica- 
tion  factor  MX  are  selected  so  that  KXM=KXMAX. 

Also,  NYE  and  MY  are  usually  selected  so  that 
KYM< <KYMAX . 

MIN, MAX  -  Real  variables  equal  to  the  minimum  and  maximum 

values  of  the  amplitude  of  the  complex  arrays  con¬ 
taining  the  aperture  fields  are  processed  by  Sub¬ 
routine  NORMH  in  preparation  for  3D  plotting  by 
Subroutine  PLT3DH . 

ROS  -  Real  variable  equal  to  the  radius  of  curvature  of 

the  outside  shape  of  the  tangent  ogive  radome . 

BOS  -  Real  variable  equal  to  the  distance  B  in  inches 

defined  in  Figure  2-1. 

FINOS  -  Real  variable  equal  to  the  fineness  ratio  of  the 

radome  as  based  on  the  outside  dimensions. 

FINIS  -  Real  variable  equal  to  the  fineness  ratio  of  the 

radome  as  based  on  the  inside  dimensions. 

The  following  variables  are  printed  when  the  receiving  patterns  are  com¬ 
puted  and  printed: 

ICUT  -  Integer  variable  which  defines  the  E-plane  (ICUT=1) 

or  H-plane  (ICUT=2)  pattern.  See  Figure  2-3. 

ICOMP  -  Integer  variable  which  defines  the  field  component 

of  the  plane  wave  incident  on  the  receiving  antenna: 

IC0MP=1  for  elevation  component;  IC0MP=2  for  azimuth 

component . 
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Real  variable  equal  to  the  sine  of  the  maximum 
angle  off  broadside  for  which  the  received  voltage 
is  computed . 

Integer  variable  (power  of  2)  equal  to  the  number 

of  points  at  which  the  receiving  pattern  is  computed. 

The  pattern  is  computed  at  NREC  points  spaced  equally 

in  k  =sinO  according  to  Ak  =2  KMAX/NREC. 
xy  xy 

Real  variable  equal  to  2 *KMAX/NREC . 

Real  variable  equal  to  sin  * (KMAX) . 

The  receiving  pattern  is  computed  at  NREC  points  and  magnified  using 
Subroutine  MAGFFT  to  256  points  equally  spaced  in  sin  0  over  the  range 
(-KMAX,  KMAX-DK) .  Three  parameters  are  printed:  angle  in  degrees,  ampli¬ 
tude  in  decibels,  and  phase  in  degrees.  Only  every  fourth  point  in  the 
256  points  is  printed.  The  receiving  patterns  are  printed  in  the  following 
order : 


E-Plane: 

£  r  - 

A 

EL 

EL 

H-Plane : 

i  , 

A  „ 

AZ 

AZ 

Subroutine  RECBS  computes  the  boresight  error  of  the  antenna  as  pro¬ 
duced  by  the  radome.  When  SUPPRS=FALSE,  the  following  parameters  are 
printed : 

Kl,  K2  -  Real  subscripted  variables  containing  the  direc¬ 

tion  cosines  (k  . ,k  , ,k  )  of  the  last  and  next  to 
xi  yi  zi 

last  true  directions  to  the  target.  One  of  these 
variables  is  equal  to  K,  the  subscripted  variable 
containing  the  direction  cosines  of  the  last  tar¬ 
get  return. 


KMAX 

NREC 


DK 

ANGMAX 


m 
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AZTM, ELTM 


MESAZ , MESEL  - 


UAZ,UEL 


SMAX 


Real  variables  equal  to  the  boresight  error  in  the 
H-plane  and  E-plane  associated  with  the  last  target 
return  (k^jk^k^) .  Expressed  in  milliradians ,  these 
errors  are  computed  according  to 


AZTM  =  Sin  1(k  //l-k  ^>*1000. 

X  y 


ELTM  =  Sin”1 (k  //l-k  2)  *1000. 
y  x 


Let  k  =  x  k  +  y  k  +  z  k  .  Then  AZTM  is  the 
Ax  Ay  A  z 

angle  between  the  z^-axis  and  the  projection  of 
k  onto  the  x  z  (azimuth)  plane.  ELTM  is  the 

A  A 

angle  between  the  z^-axis  and  the  projection  of 

k  onto  the  y  z.  (elevation)  plane. 

A  A 

Real  variables  equal  to  the  monopulse  error 
slopes  in  the  azimuth  and  elevation  channels  ex¬ 
pressed  in  units  of  volts  per  degree,  where  the 
maximum  signal  received  by  the  sum  cnannel  is  con¬ 
sidered  to  be  one  volt. 

Real  subscripted  variables  equal  to  the  received 

tracking  functions  I  {A/E}  corresponding  to  the 
mag 

target  returns  K1  and  K2  above;  e.g.,  UAZ(l)  = 

1  {A„/E  )  for  Kl. 

mag  AZ  AZ 

Real  variable  equal  to  the  maximum  amplitude  of 


tne  received  sum  channel  voltage. 


LCTR  -  Integer  variable  equal  to  the  number  of  iterations 

(target  returns)  used  by  Subroutine  RECBS  to  com¬ 
pute  boresight  error. 

Subroutine  RECBS  also  computes  and  prints  six  additional  mono pul se 

outputs  around  the  apparent  boresiqht  direction  k  .  The  directions  k  chosen 

o 

lie  in  the  plane  k  =k  and  are  spaced  one  mil li radian  apart  over  the  range 
x  y 

13  mrad  and  centered  on  the  direction  k  .  The  variables  prin'od  are  a=; 


follows : 


o 


ANG 


VRAZ , VREL 


DAZ , DEL 


SLPAZ , SLPEL  - 


Real  variable  equal  to  the  angle  in  milliradians 

between  K  and  k. 
o 

Real  variables  equal  to  I  (A/)  f  for  the  target 

mag 

return  from  direction  k  for  the  azimuth  and  ele¬ 
vation  channels,  respectively. 

Amplitude  and  phase  (degrees)  of  the  complex  volt¬ 
ages  received  on  the  A  and  A  channels,  respec- 

AZ  EL 

tively,  for  target  return  k. 

Average  values  of  the  monopulse  error  slopes 
(volts/degree)  in  the  azimuti,  and  elevation  chan¬ 
nels,  respectively,  obtained  by  a  linear  approxi¬ 
mation  of  the  tracking  functions  based  on  their 
values  at  ANG  =13  mrad.  For  example. 


SLPAZ  =  [VRAZ ( 3  mrad)  -  VRAZ(-3  mrad) ]/(. 006*57 . 3) 


The  main  program  always  prints  the  boresight  error  in  azimuth  (BSEAZ) 
and  elevation  (BSEEL) ,  and  the  values  printed  are  identical  to  AZTM  and 

.).*■ 
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ELTM  defined  above.  Main  also  computes  the  gain  of  the  antenna  in  deci¬ 


bels  with  the  radome  in  place  according  to 

GAIN  =  20.  *  ALOG10 (SMAX/VAIRM) 

For  other  than  an  "air  radome",  GAIN  is  negative  and  indicates  a  loss  in 
antenna  maximum  gain  due  to  radome  reflections  and  ohmic  (tan  6)  losses. 

The  amplitude  of  received  sum  voltage,  VAIRK,  is  always  printed  as  the  last 
item  prior  to  termination  of  the  program. 

2-4.  Comments  and  Method 

a.  Method.  The  method  of  analysis  has  been  presented  in  Section 
2.1.  Additional  details  of  analysis  are  presented  in  the  descriptions  of 
each  subroutine,  especially  Subroutine  RECM. 

b.  Supporting  Subroutines.  Twenty  eight  supporting  subroutines 

are  required  by  SIIRACP,  22  of  which  are  identical  to  those  used  by  RTFRACP 
The  purpose  of  each  one  is  briefly  described  below.  Those  subroutines 
peculiar  to  SIIRACP  and  explained  in  Chapters  3-8  are  denoted  by  asterisks 
below. 

*(1)  T RE CNF — Computes  complex  vector  aperture  electric  fields 
of  antenna  for  all  three  monopulse  channels  at  NX  x  NY 
sample  points. 

(2)  ORIENT — Computes  matrices  ROTATE  and  TRANSLate  used  for 
coordinate  transformations  by  Subroutines  POINT  and  VECTOR. 

(3)  POINT — Transforms  a  point  P(x  ,  y  ,  z  )  in  antenna  system 

AAA 

to  the  same  point  P(x  ,  y  ,  z  )  in  radome  coordinate  system, 

R  R  R 

and  vice  versa. 

(4)  VECTOR — Transforms  a  vector  from  radome  to  antenna  coordi¬ 
nate  system,  and  vice  versa. 


2' 


( 5)  X  NCPW — Comput  es  the  rectangular  •  •  1  •■ct  r  1  i-  field  component 


*(6) 

(7) 

(«) 


*(9) 

(10) 

(11) 

(12) 

(13) 

*(14) 

*(15) 


of  a  plane  wave  incident  fror..  the  direction  k  in  antenna 

A 

coordinates.  The  power  density  of  tin  t  lane  wave  m  unity. 

RECM — Computes  the  voltage  received  by  each  channel  of  tin 

antenna  for  a  plane  wave  EINC(E  ,  E ,  I.  )  incideiit  on  the 

x  y  ;• 

radome  from  the  direction  KA ( k  ,  k  k  )  . 

x  y  z 

OGIVEN--Computes  the  unit  inward  normal  to  the  tangent  ogive 
radome  surface  at  a  specified  joint. 

RXMIT- -Computes  the  transmitted  electric  fields  of  the  plane 
wave  traveling  in  direction  -k  and  incident  on  a  flat  di¬ 
electric  wall  with  unit  inner  normal  n.  The  unit  vectors 
k,  n  are  used  to  resolve  the  incident  ;  lane  wave  into  vector 
components  perpendicul ar  and  parallel  to  the  plane  of  inci¬ 
dence,  and  to  determine  the  angle  of  incidence. 

WALL-Computes  the  normal  voltage  transmission  coefficients 
of  flat  panel  model  of  the  radome  wall  as  function  of  the 
sine  of  the  incidence  angle. 

AXB — Computes  real  vector  cross  product  C  =  A  X  B. 

CAXB--Computes  the  complex  vector  cross  product  C  =  A  X  B, 
where  A  is  complex  and  B  is  real. 

RECBS — Computes  boresight  errors  of  antenna  enclosed  by  the 
radome  for  the  specified  orientation,  fineness  ratio,  etc. 

RECPTN — Computes  receiving  patterns  of  all  three  channels. 

APINT — Computes  the  fields  of  specified  planar  aperture 
fields  using  equivalent  currents. 

DIPOLES — Computes  the  fields  of  electric  and  magnetic  di¬ 
poles  located  on  a  planar  surface  as  required  by  Subroutine 
APINT. 
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*(lb)  CAXCB- -Computes  the  complex  vector  product  C  -  A  x  Is, 
where  A  and  B  are  complex. 

(17)  FAR--Computes  the  amplitude  of  the  power  pattern  from  the 

complex  plane  wave  spectra  A  (k  ,k  ) ,  A  (k  , k  )  of  an 

x  x  y  y  x  y 

antenna . 

(IB)  AMPHS- -Converts  a  compLex  num-r  from  rectangular  ti  plar 
form.  Tins  subroutine  utilizes  the  intrinsic  function 
ATAN2.  The  am[  litude  produced  is  linear  (not  decibels)  , 
and  the  phase  is  in  degrees  on  the  range  (-lHo,  18<.  )  . 

(19)  DBPV--Conver t s  a  real,  two-dimensional  array  from  linear 
to  logarithmic  values  in  decibels  or.  the  range  t  -4  IP 

(20)  NORMH--Normalizes  a  two-dimensional  real  array  to  values 
between  0  and  1. 

(21)  CNPLTH--Plots  single  dimensional  far  field  patterns  on  axi 
patterned  after  standard  pattern  recorder  paper.  CNPI.TH 
calls  Subroutine  PSI  in  addition  to  the  usual  Cal  com;  sub¬ 
routines  . 

(22)  PSI--Used  by  Subroutine  CNPLTH  to  compute  the  azimuthal 
angle  1' . 

(23)  PLT3DH--Yi e 1 ds  three-dimensional  ] lots  of  tin  data  u.  tin 
two-dimensional  real  array  FIELD.  PLT5PH  alls.  Sul  rout  in* 
PLTT,  NORMH  as  well  as  the  usual  Cal  comp  subroutines. 

(24)  PLTT-Used  by  Subroutine  PLT3DH  to  eliminate  moving  tin 
{x?n  for  hidden  lines. 

(2ri)  FFTA--Computes  the  Fast  Fourier  Transform  of  a  one¬ 
dimensional  complex  array  having  2**IEXP  elements.  Prop,  i 
operation  is  machine  dependent. 
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(26)  MAGFFT — Provides  increased  resolution  of  a  sampled  function 


using  FFT  and  Discrete  Fourier  Transform  techniques. 

(27)  JOYFFT — Provides  increased  resolution  of  selected  portions 
of  a  two-dimensional  Fourier  transform.  JOYFFT  calls  Sub¬ 
routines  FFT A  and  PWRTWO. 

(28)  PWRTWO--Used  by  Subroutine  JOYFFT  to  ensure  that  a  given 
integer  is  a  power  of  2. 

Program  Flow 

For  the  following,  refer  to  the  program  listing  in  Section  2-8  and 
line  numbers  shown  on  the  right-hand  margin  of  that  listing. 

Line  Nos.  Explanation 

1’  All  variables  bog inning  with  the  letter  K  in  the 

main  program  arc  real. 

1>  -  ,|  Declare  variable’s  and  array  dimensions.  Note  equi¬ 

valence  statements  in  Lines  21-.’>.  The  dimension  of 
IB  IT  in  Lin*  ."••*  may  1*.  computer  system  dependent . 

Not.**  it.  Line  *1  ilia'  only  twenty  fineness  ratios, 
scan  i  lain-;;,  and  scar,  angles  can  lx  accommodated. 

t.  -  Lab*-1  common  is  used  as  a  onv*ni>  nt  meant;  to  trans¬ 

mit  variat  l.*s  to  sul  r**u '  i  n<  ..  not  direct  ly  called  by 
MAIN.  Th*  1  al  **ls  ar*  g*  n*- rated  f rom  t  In  nan*  of 
!  lie  *;ul  r*  *ut  i  n<  s.  which  r — *  iv*-  t  In  van  al  1*  *,  and  *  ret 
label  i  : .  t  *  -  rmi  na  t  *  ■<!  with  t  (**.  1  *  t  t  .  r  '  *  ■  d*  :  t  • 

*.).,  HP-. ’I  •  d.-noti  •;  varialles  c.  .mm-r  t  <  •  MA 1 N  at.d 
Sul  r< >ut  i  in  Ris  M. 


nsm*  v 


40-42 


Declare  namelists  for  printing  data.  These  name- 
lists  are  no  longer  used  except  for  occasional  de- 


43-57 

61-62 

63-64 

65-67 

68 

60 


71-75 
76-80 
81-111 
112 
1  1  3 

1  14 
1  1  7 


i  i  -i-i  r< 


l.:  *-l.  l 
i  1 


bugging  purposes. 

Set  data  in  DATA  statements  as  described  above  in 
Section  2-3. 

Set  SMAX  and  VMAX  to  unity  to  prevent  division  by  zero. 
Read  and  write  TITLE  according  to  18A4  format. 

Read  input  data  using  frec-field  format. 

Compute  sine  of  the  offset  anqle  6 

OS 

Set  TABLE=FALSE  so  that  normalizing  factor  VAIRM 

can  be  computed  via  a  call  to  Subroutines  RE CM  and 

RXMIT.  In  the  latter,  TAB LE= FALSE  causes  T  ,  T.,  to 

i  h 

be  set  unity  as  in  the  case  of  no  radome. 

Write  input  data. 

Read  input  data  and  set  VAIRM  needlessly. 

Comments  explaining  input  variables. 

Set  NN=N+1“  Number  of  wall  layers  plus  one . 

Initialize  DINCH=  total  thj.ekno.ss  of  radome  wall 
in  inches. 

Read  wall  data  and  compute  total  thickness. 

Compute  DIAIN  inside  has--  d  l  am*  t  •  r  of  the  radon* 
i  n  i  nohe:; . 

>  >  'input  .  l  nd  *  .  i  of  t  h*  cent  or  . ■  i  e-m*  nt  ■ :  >;•  at  -f  i  •  i  : 

arrays  eorres:  ondi nq  to  x  - v  1  _ 

A  ‘  A 

Writ*  array  dimensional  data. 

Read  fineness  ratios,  s;  an  |  l.iii'  s,  and  i;  u  . 


u 


130-133 


Compute  wavelength  in  inches  and  centimeters. 


r 


134 


1 35-148 


140-133 

154-160 


161-168 


160-187 


Net., 


1  o  : 


Compute  3=2  7i/ X 

cm 

Compute  DAPWL=  diameter  of  antenna  aperture  in 
wavelengths. 

Convert  variables  in  inches  to  centimeters  for  input 
to  subroutines.  Some  variables  are  multiply  defined 
to  avoid  conflicts  in  labeled  common;  e.g.,  ZBOT  and 
Zl.  Note  that  DIACM  is  the  inside  diameter  of  the 
radome  in  centimeters . 

Convert  angles  from  degrees  to  radians  using  RAD=n/180. 
Compute  near  fields  of  three  channel  monopulse  antenna 
using  Subroutine  T RE CNF . 

Set  KYMAX=KXMAX,  compute  magnified  folding  wavenumbers 
KXM,  KYM,  and  print  results. 

Initialise  Cal comp  i letter,  if  required.  The  commented 
initialization  (Lines  175-185)  applies  to  the  IBM  3033 
system  at  JHll/AFL. 

Lines  18R-?(n<  ar<-  used  to  plot  tl,e  near  fields  of  the 
antenna  and/or  th,  transmitting  irincipnl  plane  power 
pa t  t erns . 

Initialize  th<  maximum  values  KMXLh,  FMXPAZ  of  the  E- 
1 1  id  M-j  l.iti"  at  t  ei-ii  ;•!  that  when  us--d  initially  as 

i  :.|  u»  t.  Lui  n  ■  it  i  t.,  PAR,  tin  r,  suit  ins  pattern  will 

ii  !i.,ra,tl  iz.e.i  with  r,  sj  •  ,  t  to  its  own  maximum  and  FMXEL 
oid  PMXLA  7.  will  t<  ...  t  •  -iu.1 1  to  these  r>  sj  ret  ive  maxima 
*.  -i) * i.!  '  FAR,  tin  resulting  patterns  will 

:•  e  r  i  i  ;  :  w:?*  t  to  FMXEL  and  FMXPAZ .  Hence 


I 


■■ 


191 

192-201 


202 


204 


205-212 


213-226 

227 

230 

233 


the  relative  gain  of  the  difference  and  sum  pat¬ 
terns  will  he  correctly  displayed  in  the  graphs. 

Iterate  for  each  of  three  monopulse  antenna  channels. 
Equate  complex  arrays  EXT,  EYT  to  the  selected  near 
field  and  compute  the  amplitude  NF  of  EXT. 

Assume  transmitting  near  fields  are  to  be  plotted 
( GRAFTR=T) . 

Call  Subroutine  PLT3DH  to  plot  the  amplitude  of  EXT. 

The  inputs  XSIZE=6.,  YSIZE=2.5,  HEIGHT=2.5  yield  a 
3D  plot  that  will  fit  on  a  8V  x  11"  report  page. 

The  inputs  NF,  NX,  NY  specify  the  real  array  to  be 
plotted  and  its  dimensions.  The  input  NMZ= . TRUE . 
directs  the  subroutine  to  normalize  NF  so  that  its 
values  be  between  0  and  1.  The  input  LDB=. FALSE 
indicates  that  the  array  NF  contains  linear  values 
rather  than  logarithmic  values  (decibels) . 

Compute  and  plot  phase  of  EXT  on  a  scale  of  -180 
degrees  to  +180  degrees.  Note  that  Line  210  ensures 
that  the  real  array  NF  contains  these  phase  values 
scaled  to  the  required  0  to  1  range. 

Repeat  amplitude  and  phase  3D  plots  for  EYT. 

Assume  GRAFSA=T  so  that  principal  plane  patterns 
are  plotted. 

If  IP=3,  go  to  Line  264  and  plot  H-plane  patterns; 
otherwise,  plot  E-plane  patterns. 

Call  Subroutine  JOYFFT  to  calculate  the  inverse  Fourier 
transform  of  the  xft-component  of  near  field  EXT  to 


33 


produce  the  plane  wave  spectrum  XEEL  from  which 
the  radiation  field  can  be  computed.  In  the  process 
of  computing  the  transform,  provide  increased  resolu¬ 
tion  from  NX  x  NY  points  to  NYE  x  NXE  points  through 
point  (NXC,NYC)  in  the  array  EXT.  In  the  k^_  direc¬ 
tion,  the  plane  wave  spectrum  is  magnified  by  MY;  it 
is  magnified  by  MX  in  the  k  direction.  The  array 

y 

FFTXY  is  a  working  array. 

234  Repeat  for  EYT  to  produce  the  plane  wave  spectrum 
YEEL  for  the  y  -comjtonent  of  field. 

235  Call  Subroutine  FAR  to  calculate  the  E-plane  elevation 
(IPWR=3)  power  pattern  FFSEL  of  the  near  field  at 
equal  samples  in  sinO  over  the  range  (-KXM,  KXM  -AK) . 
If  FMXEL<0  (and  it  is  for  IP=1) ,  normalize  FFSEL  with 
respect  to  its  own  maximum. 

237  Call  Subroutine  DBPV  and  convert  tile  power  pattern 

to  decibels  on  a  scale  of  0  to  -40  dB. 

238-241  Scale  the  values  in  FFSEL  to  the  range  of  0  to  1  for 

plotting . 

242  Call  Subroutine  CNPLTH  and  plot  the  power  pattern. 

If  KXM<1,  the  pattern  is  plotted  over  the  angular 

range  corresponding  to  sin  ^ (KXM) ;  if  KXM>1,  the 

angular  range'  is  (-90°,  90°).  Subroutine  CNPLTH 

actually  plots  conical  cuts  corresponding  to 

k  =  constant  or  k  =  constant  as  specified  by  inputs 
x  y 

KXC,  KYC.  In  the  call  here,  KXY=KYC=0  so  that  a 
principal  pattern  is  produced. 

14 


A 


243-247 

248 

249-253 


254-269 

271 

272 

273-277 

278 

279-283 

284-309 


Write  a  figure  title  for  the  plot  and  establish  a 
new  origin  for  the  next  plot. 

If  IP=2,  the  E-plane  patterns  are  finished. 

Since  JOYFFT  changes  the  input  arrays  EXT,EYT  it  is 
necessary  to  recompute  them  so  that  increased  re¬ 
solution  can  be  obtained  in  the  plane  wave  spectra  in 
the  H-plane. 

Repeat  computation  and  plotting  for  H-plane  power 
patterns. 

Iterate  the  radome  analysis  for  NFINE  fineness  ratios. 
Set  FINE  =  outside  fineness  ratio. 

Calculate  and  write  R  ,  B,  F  ,  F  as  defined  in 

OS  OS  IS 

Figure  2-1  for  the  radome  geometry. 

Compute  RDML  =  distance  from  tne  base  of  the  radome 
to  the  theoretical  tip  on  the  inside  of  the  radome. 

If  ZTOPIN<RDML,  the  radome  has  a  metal  tip,  and  a 
message  is  written  to  that  effect. 

Compute  parameters  needed  by  Subroutine  OGIVE  to 
describe  the  radome  shape.  R  and  B  are  in  centimeters 
and  apply  to  the  inside  dimensions.  AP,  the  height 
of  the  cylinder  in  centimeters,  is  not  used.  RTSQ= 
square  of  the  radius  of  the  top  disk.  RBS£>=  square 
of  the  radius  of  the  bottom  disk  (bulkhead) .  The 
other  variables,  BSQ ,  RINV,  RSQ1 ,  RP ,  and  RP2 ,  arc 
precalculated  hero  to  speed  later  computations  in 
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OGIVE. 


310 


Compute  conversion  factor  DI’MK  for  converting  milli- 


111-31*} 


315-316 

317-319 

320-322 

323-334 

333 

336-338 
3  39 


340-343 


343 


344 


radians  to  degrees. 

Initialize  the  "last"  values  of  boresight  error  in 
azimuth  (AZI.)  and  elevation  (ELL)  and  the  "last" 
value  THL  of  scan  angle.  These-  variables  are  used 
later  to  compute  boresight  error  slope  in  degrees 
per  degree  from  the  present,  and  last  values  of  bore- 
sight  error. 

Write  title  for  analysis  results. 

Write.-  parameters  of  radome  wall. 

Write  heading  for  table  of  boresight  error  and  gain 
data. 

Write  this  same  data  to  logical  unit  7  for  subsequent 
storage  as  a  disk  file,  if  desired. 

Iterate  the  radome  analysis  for  NPHI  scan  planes. 
Compute  ifi^  in  radians  as  required  by  Subroutine  ORIENT 
Iterate  the  analysis  for  NTHE  scan  angles  in  each 
scan  plane . 

Compute  6^  in  radians  as  required  by  Subroutine  ORIENT 
Call  Subroutine  ORIENT  and  compute  the  rotation  matrix 
ROTATE  and  translation  matrix  TRANSL  required  for  coor¬ 
dinate  transformations  using  Subroutines  POINT  and 
VECTOR. 

On  the  first  iteration,  TABLE  is  false  so  that  the 
maximum  amplitude  of  the  received  voltage  on  che  sum 


channel  is  computed  without  the  radome. 


345-347 


Set  the  direction  cosines  of  the  incident  plane  wave- 


348 


349-354 


355-356 


360 

362 


363 


366 


368 

369 


so  that  it  arrives  from  the  z  direction. 

A 

Call  Subroutine  INCPW  and  compute  the  rectangular 
components  PWI  of  the  incident  plane  wave  having 
polarization  specified  by  IOPT. 

Set  TSUP=T  and  TABLE=F  so  that  an  air  radome  wall  be 
used  and  so  that  printing  by  Subroutine  RXMIT  and  RECM 
will  be  suppressed. 

Call  Subroutine  RECM  and  compute  the  complex  voltages 
VR  received  on  the  sum,  difference  elevation,  and 
difference  azimuth  channels,  respectively,  correspond¬ 
ing  to  VR(I) ,  1=1,3. 

Compute  VAIRM= | VR( 1) | . 

Set  TABLE=T  so  that  on  subsequent  iterations  VAIRM 
will  not  be  recomputed,  and  so  that  the  table  of 
transmission  coefficients  will  be  utilized  when  RXMIT 
is  called. 

If  SUPPRS=F,  compute  and  print  the  E-plane  and  bi¬ 
plane  receiving  power  patterns  of  the  antenna  with 
the  radome  in  place. 

Iterate  in  J  for  E-plane  (ICUT=1)  and  H-plane  (ICUT=2) 
patterns . 

Set  the  desired  far  field  component. 

Set  the  temporary  logical  variable  TSUP=T  so  that 
printing  will  be  suppressed. 
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370-371 


372-375 


376-377 


378 


379 


380-384 

385 

386 
388 


389 


390-391 

392 


393-404 


Call  Subroutine  RECPTN  and  compute  the  complex 
received  vol tuges  on  each  of  throe  channels  at  NREC 
] joints  over  the  range  (-KMAX,  KMAX  -  DK)  . 

Increase  the  resolution  and  print  results  for  all 
three  channels.  Do  not  print  results  that  are  known 
to  be  identically  zero. 

Transfer  the  received  voltage  into  a  one-dimensional 
array  VREC. 

If  NREC>NXE,  there  is  no  need  to  increase  the  resolu¬ 
tion  . 

Call  Subroutine  MAGFFT  to  increase  the  resolution  of 
VREC  from  NREC  points  to  NXE  points.  The  result  is 
contained  in  complex  array  XYFFT  on  output. 

Compute  linear  power  pattern. 

Select  NXX=  larger  of  NXE  and  NREC. 

Write  heading  for  printed  results  from  Subroutine  NORMH . 
Call  Subroutine  NORMH  to  normalize  the  NXX  values  in 
real  array  MVREC  to  be  between  zero  and  one.  The 
input  argument  I.DR= . FALSE .  since  the  values  are  not 
in  decibels. 

Call  Subrout itu  DBPV  to  convert  the  power  pattern 
in  MVREC  to  decibels. 

Write  correct  heading  for  E-plane  or  H-plane . 

Compute  the  i ncroment  in  sinO  at  which  power  pattern 
has  been  computed  and  resolved. 

Scale  the  jxrnver  pattern  to  have1  values  between  0  and 
1.  If  SUT'PRS  "F ,  compute  the  angle  U=ANG  and  the  phase 

38 


i 


V  . 


I 


405 

406-416 


417-419 


420 


421-422 


423-425 


426-428 


of  the  pattern,  and  print  the  results  for  every 
fourth  angle. 

If  GRAFRV=T,  plot  the  receiving  power  patterns. 

Call  Subroutine  CNPLTH  and  plot  the  receiving  patterns 
in  turn.  Write  an  appropriate  figure  title  following 
each  pattern  plot.  Re-origin  the  plotter  pen  for  sub¬ 
sequent  plots.  The  result  of  Lines  330-383  is  four 

principal  plane  patterns:  E-plane  sum.  E-plane  A  , 

EL 

H-plane  sum,  H-plane  A  . 

Call  Subroutine  RECBS  and  compute  the  boresight  errors 
AZT,  ELT  in  the  azimuth  and  elevation  planes  of  the 
antenna  as  caused  by  the  radome.  On  output,  the  real 
array  KA  contains  the  direction  cosines  of  the  last 
target  return  and,  hence,  gives  the  true  direction  to 
the  target  at  the  time  that  the  tracking  functions  in 
the  azimuth  and  elevation  planes  indicated  the  elec¬ 
trical  boresight  direction. 

If  this  is  the  first  iteration  in  scan  angle,  do  not 
attempt  to  compute  boresight  error  slope. 

Compute  boresight  error  slope  (degrees/degree)  in 
azimuth  and  elevation  channels. 

Set  the  "last"  values  of  boresight  errors  and  scan 
angle  to  the  current  values  in  preparation  for  next 
iteration . 

Compute  loss  in  maximum  gain  of  the  antenna  sum 
channel  due  to  the  radome. 
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429-430  Write  results  to  logical  units  6  and  7. 

434-435  Write  maximum  amplitude  of  received  sum  voltage 

VAIRM  without  radome. 

436  Terminate  plotting  software. 

STOP 

END 

2-6.  Test  Case 

A  test  case  has  been  delivered  to  JHU/APL  under  separate  cover. 

Typical  input  data  are  shown  in  Table  2-1. 

2-7.  References 

1.  0.  K.  Huddleston  and  E.  B.  Joy,  "Development  of  Fabrication  and 
Processing  Techniques  for  Laser  Hardened  Missile  Radomes:  Radome 
El  '.ctrical  Design  Analysis",  Martin  Marietta  Purchase  Agreement 
573712,  April  1977. 

2.  G.  K.  Huddleston,  H.  L.  Bassett,  and  J.  M.  Newton,  "Parametric 
Investigation  of  Radome  Analysis  Methods",  IEEE  AP-S  Symposium 
Digest,  pp.  199-202,  May  1.978;  also,  Proc.  Fourteenth  Symposium 
on  Electromagnetic  Windows,  pp.  21-28,  June  1978. 

3.  E.  B.  Joy  and  G.  K.  Huddleston,  "Radome  Effects  on  the  Perfor¬ 
mance  of  Ground  Mappiing  Radar,"  U.S.  Army  Missile  Command,  DAAH- 
01-72-C-0598,  March  1973. 

2-8.  Program  Listing:  See  following  pages. 


C  THIS  RADOME  ANALYSIS  COMPUTER  PROGRAM,  SIIRACP,  WAS  PREPARED  FOR 
C  JOHNS  HOPKINS  APL  BY  G.K.  HUDDLESTON,  JANUARY  1980,  UNDER  THE 
C  COGNIZANCE  OF  ROBERT  C.  MALLALIEU. 

C  SUBR  TRECNF  COMPUTES  NEAR  FIELDS  SUITABLE  FOR  THIS  SURFACE 
C  INTEGRATION  (INNER  RADOME  SURFACE)  ANALYSIS  APPROACH. 

C  COMPUTED  RESULTS  ARE  ALSO  WRITTEN  TO  TAPE7  FOR  LATER  USE. 

C 

C - 

C 

C  IMPLEMENTATION  AT  APL/JHU  2/11/80  FOR  IBM  3033 

C 

C - 

C  ***  LIBRARIES  LSIIRAC  AND  MISCFFT  ARE  REQUIRED  FOR  EXECUTION  *»* 
PROGRAM  SIIRACP (INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT.TAPE7) 
IMPLICIT  REAL (K ) 

REAL  NF(4,4) ,MVREC(256)  ,KA(3) 

COMPLEX  SUMX ( 4 , 4 ) , SUM Y ( 4 , 4  ) , DELX (4,4),DELY(4,4) 

COMPLEX  DAZX(4,4) ,DAZY(4,4) .EXT (4, 4) ,EYT(4,4) 

COMPLEX  VR(16),VREC3(6,3).VREC(32) 

REAL  FFS(256,1),FFSEL(1,256) 

COMPLEX  XE (256 , 1  ) , YE (256 , 1  ) , XYFFT(5 1 2 ) , PWI ( 3 ) 

COMPLEX  XEEL(1 ,256) ,YEEL(1 ,256) 

EQUIVALENCE(XE( 1 , 1 ) ,XEEL(1 , 1 )) 

EQUIVAI,ENCE(YE ( 1  ,  1 ) ,  YEEL(  1 , 1 ) ) 

EQUI VALENCE (FFS(1 , 1) ,MVREC(1 ),FFSEL(1 , D) 

C 

LOGICAL  GR AF3D, GRAFSA .GRAFTR.GRAFRV, TABLE, SUPPRS.TSUP, SQUARE 
INTEGER  IBUF (512) 

C 

REAL  ROTATE (3. 3) ,TRANSL(3) .TITLE (18) 

REAL  FINR(20) , PHI (20) .THETA (20) 

COMMON /RECIC/DSTH , DSPHI , NTHMIN , NPHIMIN , AREA , NPOINTS , ROS , RIS , 
$ZBOTCM , ZTOPCM , BCM , RR 
COMMON /TDISKC/ZTOP , RTSQ 
COMMON /TR ACC /Z2 , Z 1 
COMMON/BDISKC/ZBOT , RBSQ 

COMMON /TRANSC/D IN (6),ER(6),TD(6),TZ , WALTOL , N , NN , D(6 ) , ZB,TK 
COMMON/OGIVC/RP, BSQ, AP, RINV, B, RSQ1 , RP2 
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21 
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38 
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C  39 

NAMELIST/GEOM/RR , RA . APIN , ZBOTIN . NX, NY , NXE , NYE , NXY , MX . MY . NXC , NYC  40 

NAMELIST/KDATA/KXMAX , KYMAX . KXM ,KYM  41 

NAMELIST/NEW/LMAX, DMRAD, IOPT , RAPMAX, VAIRM  42 

C  BOUNDARY  VALUES  NEEDED  BY  SUBR  TRACE  (INCHES,  CONVERT  TO  CM  BELOW)  4 3 

C  Z 1 =ZR  COORDINATE  OF  BOTTOM  DISK  44 

C  Z2=ZR  COORDINATE  OF  TOP  DISK  (Z1.Z2  IN  CM)  45 

C  APIN  IS  HEIGHT  OF  CYLINDER  IN  INCHES,  CONVERT  TO  CM  BELOW  46 

DATA  APIN/O . /  47 

C  ZBOTIN  IS  ZR  COORD  OF  BOTTOM  DISK  (BULKHEAD)  IN  RADOME  COORD  IN  INCHE  48 

DATA  ZBOTIN/O.OO/  4y 

C  KXMAX, KYMAX  ARE  OUTPUTS  OF  NEAR  FIELD  SUBR  50 

C  INITIALIZE  CONSTANTS  51 

DATA  RADIUS/1 EO/  52 

DATA  THETAA.PHIA.AGAM3A/0. 0,90. 0,0.0/  53 

DATA  PI/3. 1415926535898/  54 

CHHIHHIH  55 

DATA  NX , NY , NYE , NXY /4 , 4 , 1 , 5 1 2/  56 

DATA  MY/1 / , NREC/6/  57 

»»»»»•*»**»  58 

59 

READ  IN  DESCRIPTION  OF  RADOME  WALL  60 

SMAX=1 . 0  61 

VMAX=1 . 0  62 

READ(5 , 6 )TITLE  63 

WRITE(6,6)  TITLE  64 

READ(5 , * )  GRAF3D , GRAFSA.GRAFTR , GRAFRV, SUPPRS, I PENCD, SQUARE  65 

260  FORMAT(4L6)  66 

READ(5 , * )  NFINE.NPHI.NTHE.DIAOS, RAIN, RRIN.ZTOPIN, FREQ, OSANG  67 

SINOS=SIN(OSANG*PI/180. )  68 

TABLE= . FALSE.  69 

C  TABLE  IS  SET  FALSE  SO  THAT  NORMALIZING  FACTOR  CAN  BE  COMPUTED.  70 

WRITE (6 , 265 )  GRAF3D, GRAFSA.GRAFTR, GRAFRV, TABLE  71 

265  FORMAT ("  GRAF3D=" , L2, "  GRAFSA:" , L2, "  GRAFTR=" , L2, "  GRAFRV=" ,L2,  72 

$  "  TABLE=", L2)  73 

WRITE(6,270)  NFINE , NPHI , NTHE , OSANG  74 

270  FORMAT(/"  NFINE:", 15,"  NPHI:", 13."  NTHETA:" , 13, "  OSANG:  ",F5.2/)  75 

READ(5 , * )  LMAX, DMRAD, IOPT, RAPMAX, VAIRM, IPOL, ICASE, N , IPWR ,KMAX, NXE  76 
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MX-NXE/NX  "7 

IF  (MX.LT.1)  MX= 1  76 

HEAD(5  ,  *  )  DSTHIN  ,  DSPHIN  ,  NTHMIN  , NPHIMIN 

IF  (VAIRM.LE.O. )  VAIRM=1.0  60 

C  DIAOSrOUTSIDE  DIAMETER  OF  BASE  OF  TANGENT  OGIVE  RADOME  61 

:  VA IR  ^MAXIMUM  REC'D  VOLTAGE  W/O  RADOME  AT  KX-0.,KYr0.  62 

6  NFINErNO.  OF  FINENESS  RATIOS  83 

I  NPHIzNUMBER  OF  SCAN  PLANES  84 

C  NTHE=NUMBER  OF  ANGLES  IN  EACH  SCAN  PLANE  85 

C  DIAINrlNSIDE  BASE  DIAMETER  OF  RADOME  IN  INCHES  86 

C  ZTOPINrZR  COORD  (IN)  OF  TOP  DISK  (METAL  TIP)  87 

C  FREQ=FREQUENCY  IN  GHZ  88 

C  GRAF3D: .TRUE .  GIVES  3D  PLOTS  OF  INCIDENT  FIELDS  ON  APERTURE  (DELETED)  89 

C  GR AFRV= . TRUE .  GIVES  SA  PLOTS  OF  RECEIVING  PATTERNS  (AZ  &  EL)  90 

C  GR AFSA= . TRUE .  GIVES  SA  PLOTS  OF  TRANSMITTING  PATTERN  WITHOUT  RADOME  91 

C  SUPPRS: . TRUE .  SUPPRESSES  THE  PRINTING  OF  NUMEROUS  RESULTS  92 

C  RAPM AX-MAX  RADIUS  OF  ANTENNA  APERTURE  IN  INCHES.  93 

C  IOPT  SELECTS  POLARIZATION  OF  INCIDENT  PLANE  WAVE:  94 

C  =1  ELEV  (VERTICAL)  95 

C  -2  AZIMUUTH  (HORIZONTAL)  96 

C  =3  RHC  97 

C  =4  LHC  98 

C  I POL  SELECTS  POLARIZATION  OF  ANTENNA  WHEN  ICASE=1:  99 

C  =  SAME  CODE  AS  FOR  IOPT  100 

C  ICASE=1  OR  2  FOR  CIRC  APERTURE.  UNIFORM  ILLUMINATION  101 

C  =3  FOR  FLAT  PLATE  WITH  SPECIFIED  ILLUM,  VERT  POL  (CASE  III)  102 

C  N=NUMBER  OF  LAYERS  IN  RADOME  WALL  103 

C  OSANG=ANGLE  IN  DEG  IN  45  PLANE  OFF  BORESIGHT  OF  FIRST  TARGET  RETURN  104 

C  USED  BY  SUBR  RECBS  IN  GETTING  INITIAL  DATA.  105 

C  IPWR=1  FOR  POWER  IN  ELEV  COMP  OF  FAR  FIELD  PATTERN  106 

C  =2  FOR  AZIMUTH  COMP, =3  FOR  TOTAL  POWER.  107 

C  DSTHIN ^SAMPLE  DISTANCE  (IN.)  ON  RADOME  SURFACE  IN  THETA  DIRECTION  108 

C  DSPHIN=  -DITTO-  PHI  DIRECTION  109 

C  NTHMIN=MINIMUM  NUMBER  OF  SAMPLES  IN  THETA  DIRECTION  110 

C  NPHIMIN=  -DITTO-  PHI  DIRECTION  111 

NN=N+1  112 

DINCH=0 .  113 

DO  5  1=1 ,N  114 


READ(5  ,  •  )  DIN',  I  )  ,ER!  I  )  ,  TDt  I  ) 

5  DINCH-DINi I )+DINCH 
DIAINrDIA0S-DINCH*2. 

NXC  =NX/2+ 1 

N YC  =  N Y/2+ 1 

WRITE ( 6 , 4  j  NX,N1 , NXE , N YE . NXY . MX , MY 
4  FORMAT I "  NX.NY,NXE,NYE,NXY,MX,MY:",7I4) 

C  READ  FINENESS  RATIOS  FOR  THIS  RUN— BASED  ON  OUTSIDE  DIMENSIONS 
DO  1j  I=1,NFINE 
1 i  READ (5 , *  )  FI  NR ( I  ) 

C  READ  ORIENTATIONS  FOR  THU  RIJN  iDEOREESi 
DO  14  IzI.NPHI 

14  READI5 , * )  PHI(i) 

DO  15  Izl.NTHE 

15  READC5 , *  )  THETA ( I  ) 

C  COMPUTE  WAVELENGTH: 

WLIN=29.97925/(FREy»2.54  ) 

WLCMzWLIN»2.54 
BETA =2 . *PI /WLCM 
DAPWL=2 . *RAPMAX/WLIN 
C  CONVERT  TO  CENTIMETER  AND  RADIANS 
ZBOTzZBOTIN«2.54 
Z1 zZBOT 

RSQMAXz (2.59*RAPMAX)#,2 

DIACM=DIAIN*2.54 

ZTOPrZTOPIN»2.54 

ZB=ZTOP 

Z2=ZT0P 

ZTOPCM=ZTOP 

ZBOTCMzZBOT 

RA=R AIN*2 . 54 

RR=RRIN*2.54 

DSTH=DSTHIN»2.54 

DSPHI zDSPHIN*2 . 54 

RADzPI/180.0 

6  FORMAT ( 1 8A4 ) 

THETAAzTHETAA*RAD 

PHIA:PHIA*RAD 
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DO  30  lp=  1 ,  3 
DO  35  1=1, NX 
DC  35  J  =  1 ,  NY 

IF  1 TP. EQ. 1 )  EXT( I , J ) =SUMX( I , J ) 

IF  IP. EQ. 1 )  EYT ( I , J ) =SUMY ( I , J ) 

It-  (IP.E0.2)  EXT ( I ,  J )  =DELX ( I ,  J ) 

IF  (IP.EQ.2)EYT(I,J)=DELY(I,J) 

IF  (IP.EQ.3)  EXT ( I , J ) =DAZX ( I , J ) 

IF  (IP.E0.3)  EYT ( I , J ) =DAZY ( I ,  J ) 

NF( I , J ) =CABS(EXT( I , J ) ) 

35  CONTINUE 

IF  ( . NOT. GRAFTR )  GO  TO  215 
C  PLOT  3D  NEAR  FIELDS  X-COMPONENTS 

CALL  PLT3DH(6. , 2. 5 , 2. 5 , NF, NX, NY, .TRUE. , .FALSE. ) 

C  PLOT  PHASE  /LSO 
DO  i*0  1  =  1,  NX 
DO  40  J  =  1  ,  NY 
NF( I , J )=0. 

CALL  AMPHS (EXT (I , J ) , RLF , AIF ) 

NF( I , J ) = ( AIF+1 80. ) /360. 

40  CONTINUE 

CALL  PLT3DH (6. , 2. 5, 2. 5, NF, NX, NY,. FALSE . , . FALSE . ) 
C  PLOT  3D  NEAR  FIELDS  Y-COMPONENTS 
DO  45  1=1, NX 
DO  45  J= 1 , NY 
NF( I , J ) =CABS(EYT (I , J ) ) 

45  CONTINUE 

CALL  PLT3DH(6. ,2. 5,2.5, NF, NX, NY, .TRUE. , .FALSE. ) 

C  PLOT  PHASE  ALSO 
DO  50  1=1, NX 
DO  50  J=1 , NY 
NF( I , J )=0. 

CALL  AMPHS (EYT (I , J) , RLF, A IF) 

NF(I,J)=(AIF+180. )/360. 

50  CONTINUE 

CALL  PLT3DH (6 . . 2. 5 , 2. 5 , NF, NX , NY, . FALSE . , . FALSE. ) 
If  (GRAFSA)  GO  TO  215 
GO  TO  30 


215  CONTINUE  *  22g 

IF  (IP.EU.3)  GO  TO  220  230 

CATC  EL  CUT  OF  SUM  231 

NOTE  THAT  JOYFFT  CHANGES  EXT , EYT .  232 

CALL  JOYFFT ( EXT , NX, NY , MY , MX , NXC , NYC , XEEL , NYE , NXE , XYFFT , NXY , 3)  233 

CAlL  JO YFFT ( E YT , NX, NY, MY , MX, NXC , NYC , YEEL, NYE , NXE , XYFFT, NXY , 3 )  234 

CALL  FAR (FFSEL, XEEL, YEEL, NYE, NXE, FREQ, KYM.KXM, RADIUS. I PWR.FMXEL)  235 

SA  PLOTS  OF  ELEVATION  RESULTS  236 

CALL  DBPV (FFSEL, NYE , NXE , 1 )  237 

DO  216  1=1, NYE  238 

DO  216  J=1 ,NXE  239 

FFSEL ( I , J I = 1 . O+FFSEL ( I , J ) /40.  240 

216  CONTINUE  241 

CALL  CNPLTH (FFSEL , NXE , KXM , 0. ,0. )  242 

CALL  SYMBOL! .5,6.5, . 140000, 39HF IGURE  TRANSMITTING  ELEVATION  PO  243 

$WER ,0. ,39)  244 

RPWR=FLOAT(IPWR)  245 

CALL  NUMBER(999. ,999. , . 14.RPWR.0. ,0)  246 

CALL  PLOT (8. 5.0. . - 3 )  24? 

IF  (IP.EQ.2)  GO  TO  30  248 

RECOMPUTE  SUMX , SUMY  FOR  JOYFFT:  249 

CALL  TRECNF ( EXT , NX , NY , 1 , IPOL, 1 , DAPWL , DXWL , KXM AX , ICASE , SQUARE )  250 

WRITE(6,219)  I PWR  251 

219  FORMAT ( "  I  POWER  OF  PATTERN^" , 12)  25? 

CALL  TRECNF (EYT, NX, NY, 1 , IPOL, ?, DAPWL, DXWL, KXMAX, ICASE, SQUARE)  253 

220  CALL  JOYFFT(EXT, NX, NY, MX, MY, NXC, NYC, YE, NXE, NYE, XYFFT, NXY, 3)  254 

CALL  JO YFFT (EYT, NX, NY, MX, MY, NXC, NYC, YE, NXE, NYE, XYFFT, NXY, 3)  255 

CALL  FAR (FFS.XE, YE, NXE, NYE, FREQ, KXM, KYM, RADIUS, IPWR.FMXDAZ)  256 

3 A  PLOTS  OF  AZIMUTH  RESULTS  257 

CALL  DBPV(FFS, NXE, NYE, 1 )  258 

DO  10  1=1 , NXE, 1  259 

DO  10  Jrl.NYE  260 

FFS(I,J)=1.0+FFS(I,J)/40.0  261 

10  CONTINUE  262 

CALL  CNPLTH(FFS,NXE,KXM,0. ,0. )  263 

226  CALL  SYMBOL! .5,6.5, . 14000, 37HF IGURE  TRANSMITTING  AZIMUTH  POWER  264 

$ , 0. , 37  )  265 

CALL  NUMBER (999 . , 999 . , . 1 4 , RPWR , 0.  ,  0 )  266 
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CALL  PLOT (8.5,0. , -3 )  267 

30  CONTINUE  268 

205  CONTINUE  269 

C  270 

DO  100  NG= 1 , NFINE  271 

FINE=FINR(NG)  272 

C  CALCULATE  INSIDE  FINENESS  RATIO  273 

RIN=FINE*DIAOS/ (SIN (PI-2. * AT AN (2. * FINE ) ) )  2 74 

ROS=RIN*2.54  275 

BIN=RIN-DIAOS/2.  276 

FINE=SQRT ( ( RIN-DINCH) *#2-BIN**2 ) /DIAIN  277 

RDML=FINE*DIAIN+APIN  278 

IF  (ZTOPIN.LT. RDML)  WRITE(6,25)  ZTOPIN  279 

20  FORMAT ( "  TANGENT  OGIVE  PARAMETERS:  ROS(IN)="  280 

$  ,F9.5,"  BOS( IN )=" , F9 . 5 , /26X , "  FINOS=" , F5 . 3 ,  281 

$  "  FINIS=" , F8. 5 , "  RINV=" , E 12. 5 )  282 

25  FORMAT(/"  THIS  RADOME  HAS  A  TOP  DISK  AT  ZTOPIN=  ",E12.5/)  283 

C  COMPUTE  PARAMETERS  NEEDED  BY  SUBR  OGIVE  284 

R=FINE*DIACM/ (SIN (PI-2. *ATAN(2 . *FINE ) ) )  285 

RIS=R  286 

TLIS=DIAIN*FINE  287 

IF  ( ZTOPIN. GT.TLIS)  ZT0PIN=TLIS  288 

ZT0P=ZT0?IN*2 . 54  289 

ZB=ZTOP  290 

Z2=ZTOP  291 

ZTOPCM=ZTOP  292 

IF  (ZTOPCM.GT.RIS)  ZTOPCM=RIS  293 

B=R-DIACM/2.  294 

BCM=E  295 

AP=APIN*2 . 54  296 

RTSQ=R**2- (ZTOP-AF ) **2  297 

IF  (RTSQ.LT.O.)  RTSQ=0.  298 

RTSQ= (SQRT (RTSQ)-B ) **2  299 

RBSQ=R**2- (ZBOT-AP) **2  300 

IF  (RBSQ.LT.O.)  RBSQ=0.  301 

RBSQ= (SQRT(RBSQ)-B)**2  302 

BSQ=B**2  303 

RINV=1./R  304 
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RSQ1 :R**2  305 

RP=RSQ 1 -BSQ  306 

R  P2=RSQ I +BSQ  307 

WRITE (6 , 20 )  RIN, BIN. FINR(NG) .FINE,  RINV  308 

309 

DPMR= 1 80 . / ( P I* 1 000 . )  310 

AZL=0 .  311 

ELL=0 .  312 

THL=0 .  313 

TL0S=D IA0S*FINR (NG )  319 

WRITE (6, 2)  TITLE .FINR(NG) .DIAOS.TLOS, FREQ, RAIN, RRIN.DAPWL, IPOL,  316 

$ICASE , IOPT  316 

DO  8  1=1, N  317 

8  WRITE (6,7)  I,DIN(I),ER(I),TD(I)  318 

7  FORMAT(2X,I3,F13.5,F10.3,F9.9)  319 

WRITE (6,9)  320 

9  FORMAT ( //"  PHI  THETA  BSEEL  BSEAZ  SLPEL  SLPAZ  GAIN"/  321 

$  "  (DEG)  (DEG)  ( MRAD)  (MRAD)  (DEG/DEG)  (DEG/DEG)  (DB)"//)  322 

WRITE (7, 2)  TITLE, FINR (NG) , DIAOS, ZTOPIN , FREQ, RAIN, RRIN.DAPWL, IPOL,  323 

$ICASE , IOPT  324 

DO  18  1=1 ,N  325 

18  WRITE (7 , 7 )  I,DIN(I),ER(I),TD(I)  326 

WRITE (7 , 9 )  327 

2  FORMAT ( 1 H 1 , 5X , "  RESULTS  OF  RADOME  ANALYSIS  USING  INSIDE  SURFACE  IN  328 

1TEGRATI0N"/ 1 8A4/"  FINENESS  RATIO=" , F8 . 5 , 2X ,  329 

2"DIAMETER=" , F8 . 5 , "  IN,  LENGTH: ", F8. 5 , "  IN."/"  FREQUENCY:",  330 

3F8.5,"  GHZ  "/  331 

4"  RA=" , F8 . 5 , "  IN.  RR=",F8.5,"  IN.  ANTENNA  D=",F8.4,  332 

5"  WAVELENGTHS"/"  IPOL:", 12,"  ICASE:" , 12, "  IOPT:", 12//  333 

6"  LAYER  THICKNESS (IN . )  ER  TAND"/)  334 

DO  100  I  PH  1  =  1 , NPHI  335 

PHIP:PHI(IPHI)  336 

PHIR:PHIP-  337 

PHIR:PHIR»Ra.  338 

DO  100  ITHE=1 ,N*.  339 

THETAL=THETA(ITHE)  340 

THETAR: 1 80. -THETAL  341 

THETA R:THETAR*R AD  342 


CALL  ORIENT ( RA , THETAA , PHIA , RR , THETAR , PHIR , AGAM3A , ROTATE , TRANSL)  3^3 

IF  (TABLE)  GO  TO  23  344 

C  COMPUTE  NORMALIZING  FACTOR:  345 

KA( 1 )  =0 .  346 

KA(2)=0.  347 

KA ( 3 )  =  1  •  348 

CALL  INCPW(KA , PWI , IOPT )  349 

TSUP=SUPPRS  350 

TABLE=. FALSE.  351 

ZTEMP=ZTOPCM  352 

ZTOPCM=DIACM*FINE  353 

IF  ( ZTOPCM.GT. RIS)  ZTOPCM=RIS  354 

CALL  RECM(PWI,KA, NX, NY, KXMAX.KYMAX, FREQ, ROTATE, TRANSL,  355 

$SUMX , SUMY , DELX , DELY , DAZX , DAZY, VR , TABLE , TSUP , RSQMAX )  356 

C  SET  ZTOPCM  BACK  TO  THE  INPUTTED  VALUE.  357 

IF  ( ZTOPCM.LT. ZTEMP)  ZTEMP=Z',’OPCM  358 

ZTOPCM=ZTEMP  359 

VAIRM=CABS(VR(1 ) )  360 

WRITE(6, 105)  VAIRM  361 

TABLE=.TRUE.  362 

23  IF  ( .NOT.SUPPRS)  GO  TO  24  363 

GO  TO  350  364 

24  CONTINUE  365 

DO  320  J  =  1 ,2  366 

ICUT=J  367 

ICOMP=IOPT  368 

TSUP= . TRUE .  369 

CALL  RECPTN(SUMX, SUMY, DELX, DELY, DAZX, DAZY, NX. NY, ICUT, ICOMP.KMAX,  370 

$NREC,VREC3, KXMAX.KYMAX, FREQ, ROTATE, TRANSL, TABLE, TSUP, RSQMAX)  371 

DO  325  MM=1 , 3  372 

ICHAN=MM  373 

IF  ( (ICUT.EQ. 1 ) .AND. (ICHAN. EQ. 3) )  GO  TO  325  374 

IF  ( (ICUT.EQ. 2) .AND. (ICHAN. EQ. 2))  GO  TO  325  375 

DO  26  1=1, NREC  376 

26  VREC(I)=VREC3(I, ICHAN)  377 

IF  (NREC.GE.NXE)  GO  TO  31  378 

CALL  MAGFFT ( VR  EC , NREC , XYFFT , NXE )  379 

DO  305  1=1, NXE  380 


\ 


I 


305  MVREC(I)=CABS(XYFFT(I))»»2 
GO  TO  33 

31  DO  32  1=1, NREC 

32  MVREC ( I ) =CABS(VREC ( I ) ) **2 

33  NXX=MAXO(NXE , NREC ) 

WRITE (6 , 306 ) 

306  FORMAT(/"  MIN  AND  MAX  VALUES  OF  REC""G  PATTERN:  "/) 

CALL  NORMH(MVREC,NXX, 1 , .FALSE. ) 

CALL  DBPV (MVREC ,NXX, 1,1) 

IF  (J.EQ.1)  WRITE(6,308) 

IF  (J.EQ.2)  WRITE(6,309) 

DK=2.*KMAX/NXX 

IM0D=4 

IF  (NREC.GE.NXE)  IM0D=1 
DO  307  1=1 ,NXX. 1 
IF  (SUPPRS)  GO  T0307 
ANG=ASIN(-KMAX+(I-1 )*DK)*180./PI 
CALL  AMPHS(XYFFT(I) .AMP.PHS) 

IF  (NREC.GE.NXE)  CALL  AMPHS(VREC(I ) ,AMP, PHS) 

IF  ( MOD ( I , IMOD ) . EQ . 0 )  WRITE(6,310)  ANG , MVREC ( I ) , PHS 

307  MVREC ( I ) = 1 . O+MVREC ( I ) /40. 

308  FORMAT(/"  REC" "G  PATTERN,  EL  CUT,  EL  COMP  (DB) :  "/) 

309  FORMATS"  REC" "G  PATTERN,  AZ  CUT, EL  COMP  (DB)  :  "/) 

310  FORMAT (F9. 1 , 5X, F8. 3. 3X, F6. 1 ) 

IF  ( . NOT . GRAFRV )  GO  TO  325 

CALL  CNPLTH( MVREC, NXX.KMAX, 0. ,0. ) 

IF  (J.EQ.1)  CALL  SYMBOL (.5, 6. 5,. 140, 43HFIGURE  RECVG  POWER  PA 

'fcTTFRN— Ff  FV  PIflNE  0  4T) 

IF  (J.EQ.2)  CALL  SYMBOL ( .5 , 6. 5 , . 1 40, 4 1HFIGURE  RECVG  POWER  PA 

$TTERN-AZ  PLANE, 0.  , 4 1 ) 

CALL  PLOT (8. 5,0. ,-3) 

325  CONTINUE 
320  CONTINUE 
350  CONTINUE 

C  COMPUTE  BORESIGHT  ERROR 
275  CONTINUE 

CALL  RECBS(SUMX, SUMY, DELX, DELY , DAZX, DAZY, NX, NY, 

$  LMAX, NS, IOPT.VR.DMRAD, ROTATE, TRANSL, FREQ, KXMAX.KYMAX, 
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386 
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398 

399 

400 
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402 

403 

404 

405 

406 

407 
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409 
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411 
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413 

414 
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V- 


$  TABLE , SINOS , KA , AZT , ELT , RSQMAX , VMAX , SMAX , SUPPRS )  419 

IF  (ITHE.EQ. 1 )  GO  TO  300  420 

SLPAZ= (AZT-AZL) *DPMR/ (THETAL-THL )  421 

SLPEL= (ELT-ELL ) *DPMR/ (THETAL-THL )  4  22 

300  AZL=AZT  423 

ELL=ELT  424 

THL=THETAL  425 

GAINM=SMAX/VAIRM  426 

IF  (GAINM.LT. IE-2)  GAINM= 1 E-2  427 

GAINM=20. #ALOG1 O(GAINM)  428 

WRITE(6 ,11)  PHIP.THETAL, ELT, AZT.SLPEL, SLPAZ ,GAINM  429 

WRITE (7,11)  PHIP.THETAL, ELT, AZT.SLPEL, SLPAZ, GAINM  430 

11  FORMAT (1X,F5.1,F6.1,F8.2,F8.2,F9.4,F10.4,F7.1)  431 

C  GRAF3D  OPTION  HAS  BEEN  REMOVED.  432 

100  CONTINUE  433 

WRITE(6, 105)  VAIRM  434 

105  FORMAT(//"  RECEIVED  SUM  VOLTAGE  WITHOUT  RADOME=",E12.5//)  435 

IF  (GRAF 3D . OR . GRAFSA . OR . GRAFTR . OR . GRAFRV )  CALL  PLOT(0. ,0. ,999)  436 

STOP  437 

END  43& 


r.3 


I 


BLOCK  DATA 
C 

COMMON/TRANSC/DIN (6 ) .  ER (6 ) , TD( 6 ) ,  TZ , WALTOL , N , NN , D(6 ) , ZB 
C 

DATA  WALTOL.TK.TZ/O. ,0. ,0./ 

C 

END 


r>4 


* 


Chapter  3 


SUBROUTINE  RECM 


3-1.  Purpose:  To  compute  the  complex  voltages  produced  at  the  terminals 
of  the  three  channels  of  a  radome  enclosed  monopulse  antenna  by  a 
plane  wave  of  specified  polarization  and  direction  of  arrival. 

3-2.  Usage:  CALL  RECM  (EINC,  KA,  NX,  NY,  KXMAX,  FREQ,  ROTATE, 

TRANSL ,  SUMX ,  SUMY,  DELX ,  DELY,  DAZX,  DAZY ,  VREC,  TABLE, 
SUPPRS,  RSQMAX) 

COMMON/RECIC/DSTH,  DSPHI ,  NTHMIN,  NPHIMIN ,  AREA,  NPOINTS , 


ROS,  RIS,  ZBOTCM ,  ZTOPCM,  BCM ,  RR 


3-3.  Arguments 


A  complex  array  of  three  elements  containing  E  , 

E  ,  E  of  the  incident  plane  wave.  See  Subroutine 

y  z 

INCPW . 

A  real  array  of  three  elements  containing  the  di¬ 
rection  cosines  k  ,  k  ,  k  of  the  unit  vector 
xA  yA  zA 

k  which  points  from  the  antenna  origin  in  the 
A 

direction  from  whence  the  plane  wave  emanates. 


NX,  NY 


The  even  integer  number  of  sample  points  in  x 


and  y^  directions  used  to  represent  the  antenna 
aperture  fields. 

KXMAX, KYMAX  -  Real  variables  which  represent  the  normalized 

folding  wavenumbers  corresponding  to  the  sample 
distances  Ax  ,  Ay  according  to  Ax  =  A/  ( 2*KXMAX)  , 

A  A  A 


Ay  =A/(2*KYMAX) ,  where  A  is  the  free  space  wave- 

A 


length  . 


FREQ 


Frequency  in  qiualu-rtz  of  tin-  monochromatic  plane 


ROTATE, TRANS L- 


SUMX , SUMY 


wht  r- 


DEI.X ,  DELY 


DAZX, DAZY 


VRF.C 


wave1 . 

Real  matrices  of  direction  cosines  and  translation 
distances  used  to  carry  out  coordinate-  transfor¬ 
mations  of  points  and  vectors  from  antenna  to  radome 
coordinate  systems,  and  vice  versa.  See  Subroutine 
ORIENT. 

Two  dimensional  (NX  X  NY)  complex  arrays  of  the  x 

and  y  vector  components  of  the  antenna  aperture- 

fields  for  the  sum  channel  of  a  throe-channel  mono- 

pulse  antenna.  The  element  at  I=NX/2+l,  J=NY/2+l , 

corresponds  to  that  at  x  =0,  y  =0  in  the  aperture. 

A  A  1 

The  general  cor resjiondence  is  nivon  by 


<  -  x  4  ( I  -  1 )  *  .‘.x  =  ( I-MIDNX)  *Ax 

A  max  A  A 


-  v  4  (J-l)*,'v  -  ( J-MIDNY)  *Ay 

A  max  'A  A 


max  A 


x  *NX/2  and  y  -  Ay  *NY/'2  . 


max  A 


Also  .see  Sul -rout  lie  TRECNF. 

Antenna  apt  tint-  fi-  ids  for  tin-  difference  ole- 
vat  i  ■  si  hami-  ■  1  . 

Antenna  a;  i-i  ♦  mv  fit  Id-:  for  tin-  differenct'  azimuth 
chanti-  1  . 

Pom)  lex  at  ray  >  -f  !  ht  ee  e  1  <  m<  nt  winch  on  output 
e.  -iit.in-  ■  -n;  1-x  »  .  tranai  vlta*i«  of  th<-  antenna 


for  the  sum,  elevation  difference,  and 
azimuth  difference  channels,  respectively. 

TABLE  -  Logical  variable  required  by  Subroutine  RXMIT: 

if  TRUE,  a  look-up  table  is  used  to  calculate  the 
transmission  coefficients  of  the  radome  wall;  if 
FALSE,  these  coefficients  are  calculated  exactly 
for  each  angle  of  incidence  specified. 

SUPPRS  -  Logical  variable  used  to  control  the  printing  of 

results  from  Subroutine  RXMIT:  if  FALSE,  a  table 
of  power  transmission  and  reflection  coefficients 
for  equal  increments  in  the  sine  of  the  incidence 
angle  is  printed.  The  phases  of  the  complex 
voltage  transmission  and  reflection  coefficients 
of  the  radome  wall  are  also  printed. 

RSQMAX  -  Real  variable  denoting  the  maximum  radius  of  the 

2  2 

antenna  aperture  such  that  any  point  (x  +y^  ) > RSQMAX 
is  omitted  from  the  summation  procedure  used  to 
compute  the  received  voltages  VREC. 

DSTH , DSPHI  -  Real  input  variables  which  specify  the  sample  dis¬ 
tance  in  the  0  and  $  directions  on  the  radome  sur¬ 
face  ;  e . g . ,  A/3 . 

NTHMIN,  -  Integer  input  variables  which  specify  the  minimum 

NPHIMIN  acceptable  number  of  samples  N^,  in  the  two 

directions;  e.g.,  N,  =4. 

^  <j>MIN 

AREA  -  Real  output  variable  equal  to  the  surface  area  of 


the  radome  included  in  the  surface  integration. 


NPOINTS 


Integer  output  variable  equal  to  the  number  of 


sample  points  on  the  radome  surface. 

ROS,RIS  -  Real  input  variables  equal  to  the  generating  radii 

of  the  inside  and  outside  surface's  of  the  tangent 
ogive  radome  shape  (Figure  1-1) . 

ZBOTCM,  -  Real  input  variables  which  specify  the  coordinates 

ZTOPCM  of  the  bulkhead  and  opaque  tip  (if  any),  respectively 

(Figure  3  —  1 )  . 

BCM,RR  -  Real  input  variables  defined  in  Figure  3-1. 

3-4.  Comments  and  Method 

a.  Subroutines  Required:  APINT,  VECTOR,  POINT,  RXMIT,  CAXB  ,  OGIVEN, 

CAXCB . 

b.  Method:  The  voltage  V  induced  at  the  terminals  of  a  linear  antenna 

R 

by  a  "received"  electromagnetic  plane  wave  E  ,  H  is  given  by  the  Lorentz 

R  R 

reciprocity  theorem  as  [1] 

e 

VkA>  "  c  J  (^r  x  %  -  £r  x  V  *  nda  (1) 

s 

where  k  is  the  unit  vector  which  points  in  the  direction  from  whence  the 
A 

plane  wave  emanates  and  where  E^,  are  the  electromagnetic  fields  of  the 
antenna  as  ('reduced  on  the  closed  surface  S  which  surrounds  the  antenna 
when  it  is  transmitting.  The  unit  vector  n  is  the  normal  to  S  pointing 
into  the  region  not  containing  any  sources,  and  C  is  a  complex  constant. 

When  the  inside  surface  of  the  radome  is  chosen  as  (closed)  surface 
of  integration,  the  source-free  volume  is  that  inside  the  radome,  excluding 
the  space  occupied  by  the  antenna;  hence ,  n  is  equal  to  the  unit  inward 


normal  n. 

is 

elemental 


to  the  inside  radome  surface.  The  surface  can  be  divided 
areas  AA.^,  and  the  received  voltage  can  be  approximated  by 


i  rtt.o 


V  (k  )  =  C  Y  7  (E  xH  -  E  xH) 
R  A  f  L  -T  -R  -R  -T 

1  m 


n.  AA 
is  lm 


(2) 


where  the  fields  are  evaluated  at  the  same  points  P,  on  the  radome  sur- 

lm 

face.  The  elemental  areas  AA,  differ,  in  general,  from  point  to  point, 

lm 

and  must  be  included  under  the  summation. 

It  is  assumed  that  the  fields  E  ,  H  on  S  with  the  radome  in  place 
are  the  same  as  those  that  would  exist  in  the  absence  of  the  radome.  They 
are  computed  at  points  P'  from  their  specified  aperture  values  E  , 
via  the  Huygens-Fresnel  principle  as  explained  in  Chapters  5  and  6.  The 
received  fields  E  ,  H  at  P'  are  computed  by  applying  the  flat  panel  normal 

K  K 

transmission  coefficients  T  T  ,,  to  the  incident  plane  wave  E. ,  H 

nl,  n||  —l  -a 

E.  x  k  /n  at  the  point  P  on  the  outside  surface  of  the  radome  that  is  co- 
— l  A 

linear  with  the  inside  point  P'  with  respect  to  the  unit  normal  n.  . 

is 

(See  Figure  1-2) . 

The  tangent  ogive  radome  surface  is  divided  into  elemental  (trape¬ 
zoidal)  areas  by  sections  made  in  the  longitudinal  (0)  and  circumferential 

(4>)  directions.  In  both  cases,  desired  sampling  intervals  AS  ,  AS  (e.g., 

0  v 

1/3)  are  specified  as  input  data.  For  the  0  direction  of  Figure  3-1 ,  the 

number  of  samples  N  is  given  by 
0 


=  MAX 


R(0  -  0  ) 

TOP  BOT 


,  N 


AS0 


0MIN 


(3) 


where  R  is  the  generating  radius  of  the  ogive  surface 


R 


FIGURE  3-1.  RADOME  GEOMETRY  FOR  DEFINING  ELEMENTAL  SURFACE 
AREA  IN  0  DIRECTION. 


(4) 


R  =  L/  sin(-n 


Tan_1(2L/D) ) 


and  where  the  other  variables  are  defined  in  Figure  3-1.  (A  minimum 

acceptable  number  of  samples  N  „  is  also  specified) .  The  angular  limits 

GMIN 

are  given  by 


.  -1 


EOT 


=  sin  <WR) 


=  sin  1  (Z  /R) 
TOP  TOP 


(5) 


(6) 


Since  N.  is  an  integer,  the  sample  interval  AS.  must  be  recomputed  as 

O  U 


as  =  R ( e 

6  TOP 


^BOT^  /'N0 


=  RAO 


(7) 


For  iteration  in  I,  a  sample  point  at  the  center  of  an  elemental  area  on 
the  radome  surface  is  specified  by 


0  =  6.  +  AO/2  +  (1-1) *A0  ,  (8) 

and  the  corresponding  coordinate  is  given  by 

Z  -  R  sin  O  (0) 

R 

The  elemental  areas  are  formed  in  the  circumferential  (if)  direction 

as  indicated  in  Figure  3-2.  Using  AS,  as  input  data,  the  number  of  Sam¬ 
'S1 

pies  Na  in  the  <p  direct  ion  is  given  by 

v 


(.1 


Vccw 


POINT  OF  INTEREST 


l\ 


CCW  FIRST 


V 


CW  LAST 


(a)  DEFINITIONS  IN  o-DIRECTION 


(b)  TYPICAL  ELEMENTAL  VOLUME 


FIGURE  3-2  DEFINITION  OF  ELEMENTAL  SURFACE  AREA 
IN  o  DIRECTION. 


2np  I 

N  =  MAX  N 

<j>  AS  4>MIN  ( 


where  p  is  defined  in  Figure  3-1  and  is  given  by 


/  2  2 

p  =  /r  -  z  -  b 

R 


Since  N  is  an  integer,  the  sampling  distance  in  $  must  be  recomputed  as 

<t> 


AS'  =  2ap/N 


The  sample  point  at  the  center  of  an  elemental  area  is  specified  in  if  by 
(for  iteration  J) 


+  ( J— 1 )  A<f> 


where  A<ti=2Tr/N,  and  where  4>  is  a  specified  initial  point  in 

<p  o 

AA  of  a  surface  element  specified  by  (6,<j>)  is  given  by 


The  area 


AA  =  R ( AZ  -  BAG) 


AO  A0 

AZ  --  Z2  -  zl  =  R t sin ( 0  +  — )  -  sm(0  -  — 


It  is  deemed  advantageous  t  <  >  s*,*t  *  alx>v»  in  Equation  (13)  to  some 

(  j 

midpoint  of  the  illuminated  surface  <>f  the  radome .  This  is  done  by  truns 
forming  the  unit  vector  k .  to  rad-mie  coordinates;  i.e., 


k  -  x  k  +  v  k  +  z  k  x  k  t  y  k  +  >>  k 

A  A  xA  A  yA  A  7  A  R  xR  >R  yR  R  zR 


Tite  angle  :  f <>i  lows 

o 


=  cos 
o 


►  k  *  +  k  ' 
>:k  vR 


The  computations;  in  i.  proceed  first  in  the  counterclockwise  (CCW) 
direction  and  then  in  the  clockwise  (CW)  direction  as  indicated  in  Figure 
3-2 (a) .  For  the  CCW  direction,  the  coordinates  of  the  sample  point  are 
divert  by  (J=l,  tJ/2+1) 


x  =  !  cos  (i  +  (T-l)AT) 

P.  o 


y  =  p  sin  (C  +  ( J—  1 )  A(J0 

H  f 


For  the  CW  direction,  there-  results 


x  -  i  cos  (  ;  -  J.'.e' 

R  o 


Yr  =  >'  ;nn  (:<  -  J'° 


where  J  is  i  ncrem'-nte.l  from  unity  t  < .  N  4  .  (The  z  coordinate  is  given 


Kyuat  i .  •  1 1  I  1 )  . ) 


For  each  elemental  area  specified,  two  tests  may  be  performed  to 

determine  if  the  contribution  of  the  fields  on  that  surface  element  should 

be  included  in  the  received  voltage.  The  first  test  consists  of  ensuring 

that  the  sample  point  (x  ,  y  ,  z  )  lies  forward  of  the  aperture  plane  of 

R  R  R 

the  antenna;  i.e.,  that  z  >0.  The  second  test  (which  may  be  disabled 

A  - 

as  deemed  appropriate)  determines  if  the  surface  element  is  directly 
illuminated  by  the  incident  plane  wave.  The  test  is  performed  by  com¬ 
puting  the  angle  of  incidence  0  according  to 


0  5 


77-85 


86-90 

91-97 

98-104 

107-106 

107-119 

120-127 

126-137 

1 36-141 

1  4  2-132 


133-196 

137-163 

164-171 

] 72-176 


Compute  Z  ,  Z  ,  c  ,  e  and  <  nsurt  that 
1  ROS  RIS  OS  IS 

surface  clement  does  not  lie  forward  of  metal 

tip  or  aft  of  bulkhead. 

Compute  N  and  \S' . 

t  >i> 

Compute  N  ,  N  ,  and 
1  1CCW  1CW 

Compute-  AZ  and  AA. 


Iterate  in  (:  CCW  first,  CW  second. 

Compute  inside  coordinates  x  ,  y 

RIS  RIS 

Compute  unit  inward  normal  n  and  apply 
illumination  test  (disabled) . 

Convert  coordinates  of  sample  point  on  surface  to 

antenna  coordinates  to  ensure?  z  >0. 

A  - 

Compute  phase  of  incident  plane  wave  at  outside 
point  (xR0S,  yR0J. ,  ^R0S>  with  respect  to  the  an¬ 
tenna  oriqin.  Adjust  phase  of  the  specified  inci¬ 
dent  ;  l.nu-  wave  and  store  temporarily  as  H 

-RP 

Compute  antenna  coordinates  of  inside  point  in 
wave  1 enqt  h:  . 

Compute  transmitted  plane  wave  E' ,  h'  at  inside 
poi  nt . 


Use  aperture  intearation  to  compute  the  trans¬ 
mitted  fields  E  . ,  H  .of  the  antenna  at  the  inside 
-Ti  -Ti 

j-oint  for  each  channel  of  the  monojjul.se  antenna. 
Disabled  statements  pertaininq  to  surface  integra¬ 
tion  usinq  the  outside  radome  surface. 

Form  the  vector  cross  products  S  =  E  X  ll',  S.,  = 
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177-179 


180-182 


Add  contribution  to  received  voltaqe  V 

Ri 


183-184 

185-186 

187-189 

190-197 


Incement  AREA. 

Increment  NCUS  =  number  of  points  omitted. 
Increment  NPOINTS . 

If  SUPPRS=. FALSE. ,  compute  and  write  total  sur¬ 
face  area,  received  voltages,  number  of  points 
used,  and  number  of  points  omitted. 


3-6.  Test  Case:  None 
3-7.  References 
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Investigation  of  Radome  Analysis  Methods",  1978  IEEE  AP-S  Sym¬ 
posium  Digest,  pp.  199-201,  May  1978. 

3-8.  Program  Listing  (See  following  pages) 
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I 


SUBROUTINE  RECM( EINC , KA , NX, NY,KXMAX,KYMAX, FREQ, ROTATE, TRANSL,  1 

$  SUMX, SUMY, DELX , DELY , DAZX, DAZY , VREC , TABLE , SUPPRS , RSQMAX)  2 

C  SUBR  RECM  COMPUTES  THE  RECEIVED  VOLTAGE  OF  AN  ANTENNA  INSIDE  A  TANGENT 
C  OGIVE  RADOME  AS  PRODUCED  BY  A  PLANE  WAVE  INCIDENT  FROM  THE  DIRECTION  9 

C  SPECIFIED  BY  KA .  THE  INSIDE  SURFACE  OF  THE  RADOME  IS  USED  AS  THE  SURFACE 
C  OF  INTEGRATION  IN  THE  RECIPROCITY  INTEGRAL,  AND  THE  NORMAL  TRANSMISSION 
C  COEFFICIENT  IS  USED  TO  TRANSFER  THE  INCIDENT  PLANE  WAVE  FROM  THE  POINT 
C  P  ON  THE  OUTSIDE  SURFACE  TO  THE  POINT  P'  ON  THE  INSIDE  SURFACE,  WHERE  8 

C  P  AND  P'  ARE  COLINEAR  WITH  THE  NORMAL  TO  EITHER  SURFACE.  9 

C  THE  CALL  TO  THIS  SUBR  IS  IDENTICAL  TO  THE  CALL  TO  SUBR  RECM  10 

C  USED  IN  THE  RAY  TRACING  FORMULATION  HOWEVER,  ADDITIONAL  VARIABLES  11 

C  ARE  NEEDED  BY  THIS  SUBR  AND  ARE  PASSED  FROM  MAIN  PROGRAM  VIA  LABEL  12 

C  COMMON/RECIC/  AS  SHOWN  BELOW.  13 

COMPLEX  ET(3),HT(3),ERP(3),HRP(3)  19 

COMPLEX  S1(3),S2(3).U,C,EINC(3),VREC(3)  ’5 

COMPLEX  SUMX (NX, NY) , SUMY (NX, NY ) , DELX (NX, NY) , DELY (NX, NY) ,  16 

$  DAZX(NX.NY) ,DAZY(NX,NY)  17 

REAL  KXMAX , KYMAX, ROTATE ( 3 ,  3 )  . TRANSL ( 3 ) , LAMBDA , NISA ( 3 )  18 

REAL  PIR(3),NIS(3),KR(3),KA(3),PT(3).PISR(3),PO(3),PTWL(3)  19 

LOGICAL  TABLE, ATOR.RTOA, SUPPRS, INIT  20 

COMMON /R EC IC/DSTH , DSPHI , NTHMIN , NPHIMIN , AREA , NPOINTS, ROS, RIS,  21 

IZBOTCM, ZTOPCM, BCM, RR  22 

NAMELIST/ ATDR/DSTH , BETA , DKX , DKY , DKXY , THTOP , THBOT , STH , NTH  23 

DATA  ATOR/. TRUE. /.RTOA/. FALSE./  29 

DATA  PI/3.19159265/  25 

DATA  ZERO/IE-6/  26 

DATA  TUPI/6. 28318530/  27 

DATA  ETA/376.991 1  185/.NDO/0/  28 

DATA  NISA/0. ,0. ,-1 . / , PT/O. ,0. ,0./  29 

AREA=0 .  30 

B=BCM  31 

NPOINTS=0  32 

NCUS=0  33 

DKX =2 . *KXM AX/NX  39 

DKY =2 . *KYMAX/N Y  35 

NXMID=NX/2+1  36 

NYMID=N Y/2+ 1  37 

DXWL- . 5/KXMAX  38 


DYWL=.5/KYMAX  39 

IF  (NDO.CT.O)  GO  TO  4  40 

C  INITIALIZE  CONSTANTS  IN  SUBR  DIPOLES:  41 

LAMBDA=29. 97925/FREQ  42 

BETA=2.«PI/LAMBDA  43 

INIT  = . TRUE .  44 

PTWLO  )=0.  45  I 

PTWL(2)=0.  46 

PTWL(3)=2.»4.»RSQMAX/LAMBDA  47 

CALL  APINT (PTWL, SUMX, SUMY, NX, NY, NXMID, NYMID, DXWL, DYWL, ET, HT , INIT)  48 

WRITE(6,3)  PTWL(3),ET,HT  49  J 

3  FORMAK"  SUBR  DIPOLES  INITIALIZED  BY  SUBR  RECI"/  50 

$"  AT  Z=2*D**2/WL=  ",E12.5,"  ET=  ".6E12.5/  51 

$30X, "  HT=  ".6E12.5/)  52 

RTD=180./PI  53 

CALL  RXMIT(HRP,ERP,KA,NISA,PT, TABLE, SUPPRS, BETA)  54 

4  CONTINUE  55 

DKXY =DKX*DKY  56 

THTOP=ASIN (ZTOPCM/RIS)  57 

THBOT  sASIN ( ZBOTCM/ RIS )  58 

STH=RIS«(THTOP-THBOT)  59 

NTHP=STH/DSTH  60 

IF ( NTHP . GE . NTHMIN )  GO  TO  10  61 

NTH=NTHMIN  62 

GO  TO  15  63 

10  NTH=NTHP  64 

15  DSTHP=STH/NTH  65 

DTH= (THTOP-THBOT ) /NTH  66 

IF  ( .NOT. SUPPRS)  WRITE(6, ATDR)  67 

C  DETERMINE  ANGLE  PHIO  OF  CENTER  OF  ILLUMINATED  AREA  ON  RADOME:  68 

CALL  VECTOR (K A, KR.ATOR, ROTATE)  69 

RAD=KR(1 )**2+KR(2)**2  *  70 

IF  (RAD. GT. ZERO)  GO  TO  16  71 

PHIOsO.  72 

GO  TO  31  73 

16  PHIO=ACOS(KR( 1 ) /SQRT(RAD) )  74 

31  DO  32  1=1,3  75 

32  VREC(I)=(0. ,0. )  76 
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C  SELECT  CIRCLE  ON  SURFACE  OF  RADOME  AT  CONSTANT  THETA  77 

C  AND  ITERATE  IN  I  78 

TH=THB0T-DTH/2.  79 

IF  UNDO. EQ.O). AND. ( . NOT . SUPPRS) )  WRITEC6.33)  80 

33  FORMAT! 3X,"THDEG", NX, "PHIDEG", 12X,"PT" ,23X,"NIR"/)  81 

DO  20  1=1 , NTH  82 

TH=TH+DTH  83 

THD=TH*RTD  84 

SINTH=SIN(TH)  85 

PIR ( 3 ) =ROS*SINTH  86 

RHOOS=SQRT(ROS««2-PIR(3)**2)-B  87 

PISR(3)=RIS*SINTH  88 

IF  ((PISR(3).GT.ZTOPCM).OR.(PISR(3).LT.ZBOTCM))  GO  TO  20  89 

RHOIS=SQRT(RIS*»2-PISR(3)**2)-B  90 

NPHI P=T  JPI*RHOIS/DSPHI  91 

IF(N PHIP . GE . NPHIMIN )  GO  TO  40  92 

NPHI=NPHIMIN  93 

GO  TO  50  99 

40  NPHI=NPHIP  95 

C  DIVIDE  THE  INNER  SURFACE  INTO  NPHI  EQUAL  PARTS  96 

50  DSPHIP=TUPI*RHOIS/NPHI  97 

NPHI2=N PHI/2  98 

NPHICW=NPHI2  99 

NPHICCW=NPHI2  100 

IF( 2 . *NPHICCW . LT. NPHI )  GO  TO  55  101 

GO  TO  60  102 

55  NPHICCW=NPHICW+1  103 

60  DPHI=TUPI/NPHI  104 

DZ=RIS*(SIN (TH+DTH/2 . )-SIN (TH-DTH/2 . ) )  105 

DA=R IS* (DZ-B*DTH )  *DPHI  106 

1 10  DO  61  J 1=1 ,2  107 

JMAX=NPHICCW  108 

IF ( J 1 . EQ. 2 )  JMAX=NPHICW  109 

C  SELECT  A  POINT  ON  INNER  SURFACE  OF  RADOME  AT  CONSTANT  PHI  110 

C  AND  ITERATE  IN  J,  FIRST  CCW,  THEN  CLOCKWISE.  Ill 

PHI=PHIO-DPHI  112 

IF  (J1.EQ.2)  PHI=PHI0  113 

120  DO  62  J=1 , JMAX  119 


IF  (J1.EQ.2)  GO  TO  41  115 

PHI=PHI+DPHI  116 

GO  TO  42  1 17 

41  PHI=PHI-DPHI  118 

42  CONTINUE  119 

PHID=PHI*RTD  120 

CPHI=COS (PHI )  121 

SPHI=SIN(PHI)  122 

PISR ( 1 ) =RHOIS*CPHI  123 

PISR(2)=RHOIS*SPHI  124 

C  THE  POINT  OF  INTEREST  ON  INSIDE  SURFACE  HAS  RADOME  COORD  PISR (XR , YR , ZR ) . 

C  CALL  OGIVEN  TO  FIND  INNER  UNIT  NORMAL  NIS  TO  RADOME  SURFACE  126 

CALL  OGIVEN(PISR.NIS)  127 

IF  UNDO. EQ.O). AND. (. NOT. SUPPRS))  WRITE(6,56)  THD, PHID, PISR , NIS  128 

56  FORMAT(2(2X,F7.2) , 6E10. 3)  129 

C  TEST  NOW  IF  THIS  POINT  IS  ILLUMINATED  BY  PLANE  WAVE  130 

C  CUS=N IS ( 1 ) *KR ( 1 )+N IS(2 ) *KR (2 )+NIS(3 ) *KR (3 )  131 

C  IF  CUS  IS  GREATER  THAN  ZERO. AREA  IS  NOT  ILLUMINATED  132 

C  IF ( CUS. GT. 0 . )  GO  TO  59  133 

C  IF ( CUS. LT. 0. )  GO  TO  65  134 

C  GO  TO  59  135 

C  CONVERT  INSIDE  POINT  PISR (XR , YR , ZR  1  TO  ANTENNA  COORD  PT (XA , YA , ZA ) :  136 

65  CALL  PO INT (PISR . PT , RTOA . ROTATE .TRANSL)  137 

C  TEST  TO  INSURE  THAT  POINT  XR.YR.ZR  IS  ILLUMINATED  138 

C  BY  THE  ANTENNA  INSIDE  THE  RADOME  139 

C  IF  ZA>0., POINT  IS  ILLUMINATED  140 

IF  ( PT ( 3 ) ■ LT . 0 . )  GO  TO  59  141 

C  COMPUTE  PHASE  OF  INCIDENT  PLANE  WAVE  AT  OUTSIDE  POINT:  142 

PIR(1 )=RHOOS»CPHI  143 

PIR (2 ) =RHOOS*SPHI  144 

CALL  POINT(PIR , PO , RTOA , ROTATE .TRANSL)  145 

PHS=AMOD(BETA*(KA( 1 ) *PO ( 1 )+KA< 2 ) »P0(2 )+KA( 3 ) *PO( 3 ) ) , TUPI )  146 

U=CMPLX(0 . , PHS )  147 

C=CEXP (U )  148 

C  ADJUST  PHASE  OF  INCIDENT  ELECTRIC  FIELD  AT  OUTSIDE  POINT  AND  STORE  AS  HRP: 

HRP( 1 ) =EINC( 1 ) *C  150 

HRP(2)=EINC(2)*C  151 

HRP( 3 ) =E INC ( 3 ) *C  152 


C  COMPUTE  ANTENNA  FIELDS  AT  INSIDE  POINT  PISR:  153 

PTWL(1 )=PT(1 ) /LAMBDA  154 

PTWL(2 ) =PT (2 ) /LAMBDA  155 

PTWL( 3 ) =PT( 3 ) /LAMBDA  156 

C  TRANSMIT  INCIDENT  PLANE  WAVE  THRU  WALL  USING  NORMAL  XMN  COEFS:  157 

CALL  VECTOR (NIS, NISA, RTOA, ROTATE)  158 

IF  ( (NDO . EQ. 0) . AND. ( . NOT . SUPPRS) )  WRITE(6,57)  THD, PHID, PT , NISA  159 

57  FORMAT (2(2X,F7.2) , 6E 1 0. 3/ )  160 

CALL  RXMIT(HRP,ERP,KA, NISA, PT, TABLE, SUPPRS, BETA)  1 6 1 

C  COMPUTE  CORRESPONDING  MAGNETIC  FIELD*ETA :  162 

CALL  CAXB(ERP.KA.HRP)  163 

DO  58  ICH=1 ,3  164 

IF  (ICH.EQ.1)  165 

$CALL  APINT (PTWL, SUMX, SUMY, NX,NY,NXMID, NYMID, DXWL, DYWL, ET , HT, INIT)  166 

IF  (ICH.EQ.2)  167 

$CALL  APINT (PTWL. DELX.DELY, NX, NY, NXMID, NYMID, DXWL, DYWL, ET,HT, INIT)  168 

IF  (ICH.EQ.3)  169 

$CALL  APINT (PTWL, DAZX.DAZY, NX, NY, NXMID, NYMID, DXWL, DYWL, ET.HT, INIT)  170 

C  SUBR  APINT  COMPUTES  HT*ETA.  171 

C  HtiHHHHHHHOIIHHH  172 

C  THE  NEXT  TWO  STATEMENTS  ARE  FOR  OUTSIDE  SURFACE  CASE.  173 

C  CALL  POYNTIN (E , H , S)  174 

C  CALL  RXMIT(ETR , HTR , STR , NIS, PISP , TABLE , BETA , ETRP, HTRP)  175 

Q  »****»#»***»»»»«»».4»»»«***«»  175 

C  FORM  CONTRIBUTION  TO  RECEIVED  VOLTAGE  177 

CALL  CAXCB(ET,HRP,S1 )  178 

CALL  CAXCB( ERP , HT , S2 )  179 

VREC(ICH)=VREC(ICH)-((S1(1 )-S2( 1 ) ) *N IS ( 1 )+ (S 1 (2 )-S2(2 ) )#N IS(2 )+  180 

$  (SI ( 3 )-S2( 3 ) ) *N IS ( 3 )  )#DA  181 

58  CONTINUE  182 

AREA=AREA+DA  183 

GO  TO  62  184 

59  NCUS=NCUS+1  185 

62  CONTINUE  186 

NPOINTS=NPOINTS+JMAX  187 

61  CONTINUE  188 

20  CONTINUE  189 

NDO= 1  190 


I 

I 

I 


IF  (SUPPRS)  RETURN  191 

PERCNT= 1 00 . * ( 1 .-FLOAT(NCUS)/FLOAT(NPOINTS) )  192 

WRITE(6,25)  AREA , VREC , NPOINTS , PERCNT, NCUS  193 

25  FORMAT(//"  SUBR  RECI:  AREA=" , E12. 5/"  VREC=" , 6E 1 2. 5/"  NPOINTS=",  199 

415,"  PERCENT=",F5. 1 NCUS=",I6//)  195 

RETURN  196 

END  197 


Chapter  4 


v 


SUBROUTINE  TRECNF 

4-1.  Purpose:  To  compute  near-field  aperture  distributions  for  four 

types  of  three-channel  monopulse  antennas:  (1)  circular  aperture 
with  tapered  amplitude  and  uniform  phase  distributions;  (2)  flat 
plate  antenna  with  a  programmed  amplitude  distribution  and  uni¬ 
form  phase;  (3)  square  aperture  with  cos  x  amplitude  and  uniform 
phase;  (4)  single  element.  Four  polarizations  can  be  selected 
for  the  circular  and  square  apertures.  The  flat  plate  antenna  is 
vertically  (y^)  polarized  only. 

4-2.  Usage:  CALL  TRECNF  (E,  NX,  NY,  ICHAN,  IPOL,  IXY,  DAPWL, 

DXWL,  KXMAX,  ICASE,  SQUARE) 

4-3.  Arguments 

E  -  Complex  array  of  NX  by  NY  elements  which,  on 

output,  contains  the  values  of  the  specified  (IXY) 
rectangular  component  (x^  or  y  )  of  the  electric 
field  distribution  over  the  specified  (ICASE) 
antenna  aperture  having  the  specified  (IPOL)  pola¬ 
rization  for  the  sixjcified  (ICHAN)  channel  of  a 
three-channel  monopulse  antenna. 

NX, NY  -  Even  integer  number  of  points  in  a  rectangular 

array  at  which  the  aperture  distribution  is  com¬ 
puted  in  the  x^  and  y^  directions,  respectively. 
The  point  I=NX/2  +  1,  J=NY/2  +  1  corresponds  to 
x^=0,  y^=0.  For  the  single  element  case,  NX=NY=2. 
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I  CHAN 

iron 


IXY 

DAPWL 

dxwl 

KXMAX 

ICASE 


Integer  control  variable-  with  values  1,  2,  or  3 

which  selects  the  sun,  elevation  difference,  or 

azimuth  di  fforet.ee  channel,  respectively. 

Integer  control  variable  which  selects  the  antenna 

polarization  as  follows; 

I  -  Vertical  (y  )  ;  olarization 

2.  -  Horizontal  (x  )  " 

A 

1  -  Right-hand  circular  " 

4  -  Left-hand  circular  11 

Integer  control  variable  having  values  1  or  2  to 

select  th-  x  or  v  component  of  aperture  electric 
A  '  A 

field. 

Diametei ,  iii  wavelengths,  of  the  antenna  aperture. 
Spacing,  in  wavelengths,  between  samples  in  aper¬ 
ture  in  x  and  v  directions  (output). 

A  '  A 

Maximum  value  of  normalized  wavenumber  correspond¬ 
ing  to  KXMAX  =  l./(2.*DXWL)  (output). 

Integer  control  variable  having  values  1  or  2  to 
specify  a  circular  aperture  antenna  with  uniform 
amplitude  and  phase.  If  ICASE-.),  a  flat  plate 
antenna  having  a  programmed  am) litudo  distribution 
(see  Table  with  vertical  polarization  is 

selected . 

Logical  input  von. it  !•-;  if  TREE,  square  aperture 
is  US'-d. 


7  b 


SQUARE 


4-4. 


Comments  and  Method 


a.  The  integers  NX, NY  must  each  be  equal  to  each  other  and  even; 
e.g.,  NX=NY=l(i .  In  addition,  when  ICASE=3  (flat  pilate  antenna),  NX  and  NY 
must  equal  16.  If  NX=NY=2,  the  fields  of  a  single  element  at  x^=yA=0  are 
specified.  If  NX=NY=32,  only  the  central  15  x  15  elements  are  non-zero. 

b.  The  actual  shape  of  the  circular  aperture,  as  approximated  by  a 
rectang  a-  array  of  sample  p>oints,  is  shown  in  Figure  4-1  for  the  case  of 
NX=NY=16.  Row  1  and  Column  1  of  the  array  contain  null  elements.  The  ele¬ 
ments  inside  and  on  the  boundary  of  the  aperture  may  contain  non-zero  values 
as  shown  in  Table  4-1  for  the  various  cases  when  ICHAN=1  (sum  channel). 

Not  that  specification  of  in  Figure  4-1  determines  the  sample  spacings 

according  to 


D  cosa 

AP 

(N  -2) 
x 


D  cosa 

_AP  - 

(N  -2) 

y 


(1) 


where  a  =  Tan  '*'(2/7). 

The  aperture  distributions  for  three  monopulse  channels  are  formed 

by  phasing  the  elements  in  the  four  quadrants  of  the  aperture  appropriately. 

The  sum  channel  distribution  is  formed  by  assigning  equal  phases  to  all 

elements.  The  azimuth  difference  channel  is  formed  by  multiplying  all 

elements  in  Quadrants  II  and  III  of  the  sum  distribution  by  minus  one  and 

by  zeroing  all  elements  along  x  =0.  For  the  elevation  difference  channel, 

A 

Quadrants  III  and  IV  are  negated,  and  all  elements  along  the  line  y  =0  are 
made  zero  for  symmetry  reasons. 

The  ['basing  chosen  models  a  tracking  antenna  and  provides  outputs 
in  two  orthogonal  channels  from  which  the  direction  of  arrival  of  a  target 
return  can  be  mathematical ly  determined.  Let  k  be  a  unit  vector  which 
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I 

] 


OF  CIRCI 
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-1.  APPROXIMATION 
BY  RECTANGUU 


V 


points  from  the  antenna  origin  toward  the  dii-ection  from  whence  the  plane 
wave  (target  return)  emanates;  i.e.. 


x  k  +  y„  k  +  z  k 
Ax  Ay  A  z 


(2) 


Define  the  tracking  functions  for  this  plane  wave  as 


f .  (k  ,k  ) 
l  x  y 


A.  (k  ,k  ) 
i  x  y 

I(k  , k  ) 
x  y 


(3) 


where  A^  represents  the  output  of  the  elevation  (t)  or  azimuth  (a)  differ¬ 
ence  channel  and  Z  represents  the  sum  channel  output.  Then  for  small  k^_>0, 

the  phase  of  f  is  +  tt/2  ;  for  small  k  <0,  the  phase  of  f  is  —  tt/ 2 .  Simi- 
a  x  a 

larly,  for  small  k  >0,  arg  (f  )  =tt/2;  for  small  k  <0,  arg  (f  )  =  -tt/2. 

y  e  y  e 

Hence,  the  change  in  phase  by  tt  in  either  channel  represents  the  boresight 
direction  of  the  antenna,  and  tracking  is  done  using  the  imaginary  parts 
of  the  tracking  functions  rather  than  their  real  parts. 

c.  The  shape  and  sampling  grid  used  to  model  the  flat  plate  antenna 
are  shown  in  Figure  4-2.  In  Subroutine  TRECNF,  the  integers  NX  and  NY  must 
both  equal  16,  and  only  linear  polarization  (y^)  is  applicable  to  the  flat 
plate  antenna  (ICASE=3).  The  phasing  of  the  four  quadrants  is  done  as 
described  above  to  model  the  three  monopulse  channels  so  that  tracking  can 
be  simulated.  Note  that  specification  of  D  determines  the  sample  spacing 
according  to 


Ax 


Ayfl  = 


D  cosa. 

AI> 


<T-  2) 


wlie  re 


Tan 


(4/6)  . 
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(4) 


Table  4-1.  Values  of  Non-2ero  Elements  in 
(ICHAN=1,  ICASE=1  or  2) 


I  POL 
1 
1 

2 

2 

3 

3 

4 


IXY 

1 

2 

1 

2 

1 

2 

1 


Value 
<0  +  j0) 
(1  +  j0) 
(1  +  j0) 
(0  +  jO) 
(0  +  jl) 
(1  +  jO) 
(0  -  jl) 


Circular  Aperture 

Polarization  Type 

Vertical 

>1 

Horizontal 

fi 

RHC 

II 

LHC 


4 


2 


(1  +  jO) 


The  phase  of  each  sample  point  in  Figure  4-2  for  the  sum  channel 


is  made  equal,  but  the  amplitudes  are  tapered  in  the  and  directions 
as  shown  in  Table  4-2.  The  amplitude  distribution  is  separable  and  sym¬ 
metrical  so  that 


E 

yA 


(VV  = 


g(XA)h(V 


E  J- 

yA 


vv 


E  (x  , -y  ) 
yA  A  A 


(5) 


It  is  noted  that  samples  10,  12,  14,  and  16  are  actually  specified  in 
Lae  program,  and  samples  9,  11,  13,  and  IS  are  obtained  from  them  by 
averaging . 

d.  The  square  aperture  is  formed  by  setting  to  zero  How  1  and 

Column  1  of  the  array  of  Figure  4-1  for  symmetry  reasons.  The  values 

of  field  at  the  other  points  in  the  aperture  are  computed  to  yield  a  cos  x 

amplitude  taper  in  the  x  direction  and  a  uniform  amplitude  in  y  ;  i.e., 

A  A 


E(x,y) 


cos 


TTX 


2  x 

max 


(5) 


where  x  corresponds 
max 

3-9.  Program  Flow 
Iji ric  Nos. 

1'? 


20-22 


23 


to  the  sample  at  I=NX. 


Comments 

Assign  complex  values  to  CFAC  to  use  in  generating 
vertical,  horizontal,  RHC,  and  LHC  polarization 
according  to  IPOL. 

Compute  the  angle  a  and  the  upper  bound  R  of 

max 

the  radius  of  the  circular  aperture. 

Ensure  that  IPOL  has  correct  values  of  1,  2,  3, 


Table  4-2.  Symmetrical  Amplitude  Distribution  for  Flat  Plate  Antenna 


Sample  No. 

X 

A 

Amplitude 

^A 

Amplitude 

9 

0 

1.0280 

0 

1.0280 

10 

Ax 

1.0280 

Ay 

1.0280 

11 

2  Ax 

.9120 

2  Ay 

.9170 

12 

3Ax 

.7959 

3Ay 

.8060 

13 

4Ax 

.6077 

4  Ay 

.6155 

14 

5Ax 

.4194 

5  Ay 

.4250 

15 

6  Ax 

.2097 

6Ay 

.2125 

16 

7  Ax 

0.0 

7  Ay 

0.0 
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If  NXj^NY  and  SQUARE=FALSE ,  write  error  message 


26-21) 

30-31 

32 

33 


34-56 


57-60 


61-63 


64-65 


66-67 


6  H 


60-72 


and  stop  the  program . 

Compute  indices  of  midpoint  17. 

Ensure  that  IXY=1  or  2. 

If  NX  and  NY  are  not  even,  stop  the  program. 

Test  value  of  ICASE:  if  ICASE=3  generate  fields 
of  flat  plate  antenna  (Lines  64-105);  otherwise, 
generate  fields  of  circular  or  square  aperture 
(Lines  34-60) . 

Assign  complex  field  value  to  each  sample  point 
(x  ,y  ,0)  in  the  circular  aperture  according  to 

A  A 

the  values  shown  in  Table  4-1.  If  /x  ^  +  y  2  > 

A  A 

R  ,  make  the  field  value  zero.  Multiply  the 
max 

non-zero  elements  by  CFAC(IPOL)  to  generate  the 
correct  polarization.  For  the  square  aperture, 
zero  Column  1  and  Row  1,  and  insert  cos  x  taper 
(Line  37) . 

Compute  sample  spacing  Ax^/A  and  go  to  statement 
60. 

Error  message  and  STOP. 

Flat  plate  antenna —  if  NX^16,  write  error  message 
and  STOP  (Lines  J 31-133). 

Compute  sample  spacing  Ax  /A. 

Ensure  NX=NY 

Zero  all  elements  in  tile  aperture.  If  IXY=1 
( x  -component) ,  to  to  statement  60. 
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73-80 


Assign  tapered  amplitude  values  to  eight  "even 


81-89 

90-93 

94-97 

98-105 

106 

107-111 

112-120 

121-129 

130-134 


elements  in  Quadrant  III. 

Compute  amplitude  values  for  the  "odd"  elements 
in  Quadrant  III. 

Compute  amplitude  values  for  elements  3-9  along 

y  =0  line  and  along  x  =0  line. 

A  A 

Generate  symmetrical  amplitude  values  in  Quadrant 
IV. 

Generate  symmetrical  amplitude  values  in  Quadrants 
I  and  II . 

Compute  k  and  d  /A. 

xmax  x 

Test  to  determine  if  the  sum  channel  data  generated 
should  be  phased  to  produce  the  aperture  distribu¬ 
tion  for  a  specified  difference  channel  (ICHAN) . 
Form  aperture  distribution  for  difference  elevation 
channel  by  zeroing  all  elements  along  yA=0  and 
negating  all  elements  for  y^cO.  RETURN. 

Form  aperture  distribution  for  difference  azimuth 
channel  by  zeroing  all  elements  along  x^=0  and 
negating  all  elements  x^<0.  RETURN. 

Error  messag  for  ICASE=3  and  NX^16 . 

END 


4-6.  Test  Case:  None. 
4-7.  References 


1.  D.  R.  Rhodes,  Introduction  to  Monopulse,  McGraw  Hill,  New  York, 
1959. 


4-8.  Program  Listing:  See  following  pages. 
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SUBROUTINE  TRECNF ( E , NX , NY , ICHAN , IPOL , IXY , DAPWL , DXWL , KXMAX, ICASE,  1 

&SQUARE)  2 

C  ***  MODIFIED  JAN  80  FOR  SQUARE  APERTURE  AND  FOR  SINGLE  ELEMENT***  3 

C  SUBR  TRECNF  COMPUTES  ELECTRIC  FIELDCOMPONENTS  OVER  A  CIRCULAR  APERTURE  4 

C  OF  RADIUS  RMAX= (NX/2-1 )/COS(ATAN(2/7 ) )  AND  RETURNS  SAME  IN  E (NX, NY ) .  5 

C  NX  MUST  EQUAL  NY  AND  MUST  BE  EVEN.  6 

C  ICHANrl  FOR  SUM  CHANNEL  IP0L=1  FOR  VERT-Y  POL.  IXY  =  1  FOR  X-COMP.  7 

C  =2  FOR  ELEV  DIFF  =2  FOR  HORIZ-X  POL  =2  FOR  Y-COMP.  8 

C  =3  FOR  AZ  DIFF  =3  FOR  RHC  POL  9 

C  =  =4  FOR  LHC  POL  10 

C  DAPWL=D IAMETER  OF  APERTURE  IN  WAVELENGTHS  (INPUT)  11 

C  DXWL=SAMPLE  SPACING  IN  APERTURE  (OUTPUT)  12 

C  KXMAXrMAXIMUM  WAVENUMBER  (OUTPUT)  13 

C  ICASE=1  OR  2  FOR  UNIFORM,  CIRCULAR  APERTURE  (ADA  M.’S  CASE  I  AND  II)  14 

C  =3  FOR  FLAT-PLATE  ANTENNA,  VERTICAL  POL  (CASE  III).  15 

COMPLEX  E (NX , NY )  , CFAC( 4  )  16 

REAL  KXMAX  1? 

LOGICAL  SQUARE  18 

DATA  CFAC/(1.,0.),(1.,0.),(0.,+1.),(0.,-1.)/  19 

ANG=ATAN(2./7.)  20 

IF  (ICASE. EQ. 3)  ANG=ATAN(4./6. )  21 

RMAX=(NX/2-1 )/COS(ANG)+.001  22 

IF  (IPOL.GT.4)  IPOL=4  23 

IF  (IPOL.LT.1)  IPOL=1  24 

IF  ((.NOT. SQUARE). AND. (NX. NE. NY))  GO  TO  15  25 

NXMM7  =NX/2+1 -7  26 

NXMP7  =NX/2+1 +7  27 

NYMM7=NY/2+1-7  28 

NYMP7=NY/2+1+7  29 

C  FOR  NX,NY=32,  ONLY  THE  CENTRAL  15  X  15  ELEMENTS  ARE  NONZERO.  30 

IF  ( (IXY.LT. 1 ) .OR. (IXY.GT.2))  IXY=2  31 

IF  ( MOD(NX , 2 ) . NE. 0)  GO  TO  15  32 

IF  (ICASE. EQ. 3)  GO  TO  25  33 

TUXMX=FLOAT (NX)  34 

DO  10  1=1 , NX  35 

X=FL0AT(-(NX/2)+I-1 )  36 

COSX=COS(3. 14159265*X/TUXMX)  37 

DO  10  J=1 , NY  38 


IF  (U.EQ.I).OR.(J.EQ.I))  GO  TO  9  39 

IF  (NX.EQ.16)  GO  TO  1  40 

IF  ( (I . LT. NXMM7 ) .OR. ( I . GT. NXMP7 ) .OR. ( J . LT. NYMM7 ) .OR. ( J . GT. NYMP7 ) )  41 

$GO  TO  9  42 

1  IF(SQUARE)  GO  TO  8  43 

Y=FLOAT(-(NY/2)+J-1 )  44 

R=SQRT (X**2+Y**2)  45 

IF  (R.GT.RMAX)  GO  TO  9  46 

8  IF  ( (IPOL.EQ. 1 ) .AND. (IXY.EQ. 1 ) )  GO  TO  9  47 

IF  ((IPOL.EQ. 2). AND. (IXY.EQ. 2))  GO  TO  9  48 

C  IF  RHC ,  EY  =  ( 1 , 0) ,  EX= (0 , 1 )  I.E.,  EX  LEADS  EY  BY  90  DEG.  49 

C  IF  LHC , EY= (1 , 0) ,  EX= (0 , -1 )  I.E.,  EX  LAGS  EY  BY  90  DEG.  50 

E(I, J)=CMPLX(C0SX,0. )  51 

IF  ( ( I  POL . LT . 3 ) -  OR . ( I XY . EQ . 2 ) )  GO  TO  10  52 

E(I, J)=E(I, J)*CFAC(IPOL)  53 

GO  TO  10  54 

9  E (I , J ) = (0 . ,0. )  55 

10  CONTINUE  56 

IF  (NX. EQ. 2)  GO  TO  56  57 

DXWL= (DAPWL/2 . ) *C0S( ANG ) /(NX/2-1 )  58 

IF(SQUARE)  DXWL= (DAPWL/SQRT (2 . ) )/(NX-2)  59 

GO  TO  60  60 

15  WRITE (6,20)  61 

20  FORMAT(//"  NX. NE. NY  OR  NX  NOT  EVENIN  SUBR  TRECNF"//)  62 

STOP  63 

C  THE  FOLLOWING  IS  FOR  ADA  M.'S  ddE  III  (ICASE=2):  64 

25  IF  (NX. NE.  16)  GO  TO  90  65 

DXWL= (DAPWL/2 . ) *COS( ANG ) /( NX/2-2 )  66 

IF(SQUARE)  DXWL=(DAPWL/SQRT(2. ) )/(NX-2)  67 

NY=NX  68 

DO  26  1=1, NX  69 

DO  26  J= 1 , NY  70 

26  E(I, J)=(0. ,0. )  71 

IF  (IXY.EQ. 1 )  GO  TO  60  72 

E(6,4)=( .2824,0. )  73 

E(8.4)=(. 4250,0.)  74 

E (4 , 6 )= ( .2888,0. )  75 

E(6,6)=(. 5218.0.)  76 


E(8,6)=( .8060,0. )  77 

E(4,8)=(. 4194,0. )  78 

E(6,8)=(. 7959,0.)  79 

E (8 , 8 ) = ( 1 .028,0. )  80 

DO  30  J =4, 8, 2  81 

DO  30  1=3.8, 1  82 

IF  ( (MOD( J ,2).EQ.0).AND.( MOD ( I , 2) . EQ. 0 ) )  GOTO  30  83 

E(I,J)=(E(I— 1 , J )+E (1+1 ,J))/2.  84 

30  CONTINUE  85 

DO  35  1  =  3.8, 1  86 

DO  35  J=3 , 8, 2  87 

E(I,J)=(E(I,J-1 )+E(I,J+1 ))/2.  88 

35  CONTINUE  89 

DO  40  1=3,9  90 

40  E(I, J  )  =E (1 . 8 )  91 

DO  45  J=3 , 9  92 

45  E(9,J)=E(8,J)  93 

DO  50  J=3 , 9  94 

DO  50  1=1,6  95 

E(9+I,J)=E(9-I,J)  96 

50  CONTINUE  97 

DO  55  1=3,15  98 

DO  55  J= 1.6  99 

E (1 , 9+J )=E ( 1 , 9-J )  100 

55  CONTINUE  101 

GO  TO  60  102 

56  DXWL=DAPWL/SQRT(2. )  103 

KXMAX= . 5/DXWL  104 

RETURN  105 

60  KXMAX= 1 . / (2 . *DXWL )  106 

IF  (ICHAN.EQ.1)  RETURN  107 

IF  ( ( IXY . EQ. 1 ) . AND. ( I CASE. EQ. 3 ) )  RETURN  108 

IF  ( ( IXY . EQ. 1 ) . AND. ( I  POL. EQ. 1 ) )  RETURN  109 

IF  ( ( IXY . EQ. 2) . AND. ( IPOL. EQ. 2) )  RETURN  110 

IF  (ICHAN.EQ.3)  GO  TO  75  111 

C  LOAD  ELEVATION  DIFFERENCE  CHANNEL:  112 

J  =N Y/2+ 1  113 

DO  65  1=1 .NX  114 


65  E(I, J ) = (0 . ,0. )  1 15 

JMAX=N Y/2  116 

DO  70  J  =  1 , JMAX  117 

DO  70  1=1, NX  118 

70  E(I, J)=-E(I,J)  119 

RETURN  120 

C  LOAD  AZIMUTH  DIFFERENCE  CHANNEL:  121 

75  I =NX/2+1  122 

DO  80  J=1,NY  123 

80  E (I , J )  =  (0 .  ,0.)  124 

IMAX=NX/2  125 

DO  85  1=1 , IMAX  126 

DO  85  J=1,NY  127 

tf5  E (I , J ) =-E (I , J )  128 

RETURN  129 

C  DAPWL=5 . 047  FOR  ADA  M.'S  CASE  III  130 

90  WRITE (6 , 95 )  131 

95  FORMAT ( //"**#ERROR  EXIT!  NX  NOT  EQUAL  TO  16  IN  SUBR  TRECNF***"/ /)  132 

STOP  133 

END  134 


Chapter  5 


SUBROUTINE  APINT 


5-1.  Purpose:  To  compute  the  electromagnetic  fields  E,  H  of  a  rectan¬ 
gular  aperture  in  the  z=0  plane  at  a  point  P(x,y,z>0),  where  the 

amplitude  and  phase  of  the  aperture  electric  fields  E  ,  E  art 

xap  yap 

specified  at  N  by  N  discrete  points  spaced  d  /A  and  d  /A  apart. 

X  y  x  y 

The  aperture  magnetic  fields  H  ,  H  are  derived  from  E  via 

xap  yap  — ap 

the  geometrical  optics  approximation. 

5-2.  Usage:  CALL  APINT  (PFWL,  EX,  EY,  NX,  NY,  MIDX,  MIDY, 


DXWL,  DYWL,  E,  H,  INIT) 


5-3.  Arguments 


Real  input  array  of  three  elements  which  specifies 
the  Cartesian  coordinates  in  wavelengths  of  the 
point  P(x/\,  y/A,  z/A)  at  where  the  fields  are  to  b. 
computed;  i.o.,  PFWL(l)  =  x/A,  etc. 


EX,EY 


Complex  in; ut  arrays  of  NX  by  NY  elements  each 


which  specify  the  aperture  electric  field. 


NX, NY 


Integer  input  variables  equal  to  the  number  of 


MIDX, 

MIDY 


DXWL, 

DYWL 


sample  points  in  the  aperture  in  the  x  and  v 
directions,  respectively.  NX  and  NY  must  he  even. 
Integer  input  variables  equal  to  the  indices  in  the 
arrays  EX,  EY  corresponding  to  x=y=<A;  i.o.,  MIDX- 


NX/2+1 ,  MlDY=NY/2+l. 


Real  input  variables  equal  to  the  sample  spacings 
in  wavelengths  in  the  and  y  direct  i  ns,  n  spre- 


t ively . 


IHECEbl 


Bt.w'.K-eOr 


Com}  lex  out  ',  ut  array;,  of  three  elements.  each  equal 
t-o  the  r«  ;  tunqular  vector  comi  on. -lit;;  of  the  elec¬ 
tric  and  magnetic  fields  at  1;  i.e.,  E ( 1 ) =E  ,  etc. 

x 

Ioyica!  i  i»f  ut  variable  which  controls  i  i.i  t.ial  iza- 
t  i  on  *>f  Subroutine  S';  I  POLES  . 
d 

y ,z)  due  to  thi'  ;imi  led  aj.ort.urt  fields  arc  con- 
individual  field:;  of  equivalent  electric:  cr.d  r,ac- 
'  all  samt  le  ;  o.  tit  as  exs  lamed  in  Section  i-4. 

Declare  van  .tides,  initialise  constants  . 
con:  ut.  in.i  tial  source  joint:,  minus  d  />. .  Set 
;:-.:oori]in.Jti  .f  source  joints  to  cere. 

Initialize  summation:  of  the  fields  F , H . 

Com;  ut-.  first  sour:’,  i  oint  V  (x,v,0)  . 

x 

,  .  ■  *.  rr.  in 

e'er’.:  ir .  ■  1 ..  t  r  i  ,i  ,  d  ai.n  maoi.et  ic  J  ,  J  currcn 

x  y  '  >:  y 

.-■c.j, i  :i  :  !  .  J  c  x  H  .me.  J"'  E  x  . 

"‘Vi  I  ‘  — cl  ? 

'oil  Sui  rout  it.-  I  r  POLKS  1 1  >  ’  'To :  u  t  •  tdi.  f  i  •  Ids  of 

el.c'ri..  and  pmci.  t  >. ..  d  i  ;  l.  s  J*  ,  j‘  local,  d 

it  tile  q  -i  ■(.’  i  i  i  .  .  i  seu  r  •  :  i  ;  :  •  . 

Add  rdr.t  ri  but  i-.  >n  ;  ■  . .  •  •  mi.ir  .  -to:  ”  i .  :  t  te 

!  field  I  .  . 


put  t  ■  ; 


11.’’ 


UNCLASSIFIED 


6E0R6IA  INST  OF  TECH  ATLANTA  F/G  17/9 

PARAMETRIC  INVESTIGATION  OF  RAOOmE  ANALYSIS  METHODS.  VOLUME  III—ETC(U> 
FEB  81  G  K  HUOOLESTON,  H  L  BASSETT  AF0SR-77-3969 

AFOSR-TR-8I-0461  NL 


SUBROUTINE  APINT ( PFWL , EX , EY , NX , NY , MIDX , MIDY , DXWL , DYWL , E , H , INIT )  1 

C  SUBR  APINT  COMPUTES  FRESNEL  FIELDS  OF  RECTANGULAR  APERTURE  WITH  2 

C  APERTURE  FIELDS  GIVEN  BY  EX,EY(  H  FIELDS  ARE  DERIVED  USING  G.O.  APPROX.) 

C  FIELDS  E,H  ARE  COMPUTED  AT  THE  POINT  PFWL.  4 

COMPLEX  E(3),H(3).EX(NX,NY),EY(NX,NY),JE(2),JM(2),ES(3),HS(3)  5 

C  JE.JM  ARE  ELECTRIC  AND  MAGNETIC  SURFACE  CURRENT  DENSITIES  FOUND  FROM  6 

C  EAPXZHAT  AND  ZHAT  X  HAP.  7 

LOGICAL  INIT  8 

REAL  PSWL(3),PFWLC3)  9 

DATA  ETA/376.9911185/  10 

C  NX, NY  MUST  BE  EVEN  SO  THAT  OMITTING  ROW  1  AND  COL  1  YIELDS  SYM  APERTURE 
C  INIT=.TRUE.  TO  INITIALIZE  CONSTANTS  IN  SUBR  DIPOLES  12 

PSWL(1 )=(1-MIDX)*DXWL  13 

PSWL(3)=0.  14 

DO  1  L=1 , 3  15 

E(L)=(0.,0.)  16 

H(L)=(0.,0.)  17 

1  CONTINUE  18 

DO  10  1=2, NX  19 

PSWL ( 1 ) =PSWL ( 1 ) +DXWL  20 

PSWL(2)=(1-MIDY)*DYWL  21 

DO  10  J  =2 , NY  22 

PSWL (2)=PSWL(2) +DYWL  23 

JE(1 )=-EX(I, J)/ETA  24 

JM( 1 )=EY(I , J)  25 

JE (2 ) =-EY ( I , J ) /ETA  26 

JM(2)=-EX(I, J)  27 

CALL  DIPOLES( JE.JM, PSWL, PFWL, DXWL, DYWL.ES.HS, INIT)  25 

DO  5  L=1 , 3  29 

E(L)=E(L)+ES(L)  30 

H(L)=H(L)+HS(L)  31 

5  CONTINUE  32 

10  CONTINUE  33 

RETURN  3^ 

END  35 
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Chapter  6 

SUBROUTINE  DIPOLES 

6-1.  Purpose:  To  compute  the  electromagnetic  fields  E  =  x  +  y  E^  + 

z  E  and  H  =  x  H  +yH  +zH  at  point  P,(x/A,  y/A,  z/A)  as  pro- 
z  x  y  z  i 

e  m 

duced  by  electric  J  =  z  x  H  and  magnetic  J_  =  E^  x  z  surface  cur¬ 
rents  flowing  on  the  planar  rectangular  surface  of  dimensions  Ax/A, 
Ay/A  located  at  source  point  Pg  ,  y'/A,  z'/A)  and  oriented 

in  the  z=z*  plane.  All  dimensions  are  in  wavelengths. 

6-2.  Usage:  CALL  DIPOLES  (JE,  JM,  PSWL,  PFWL,  DXWL,  DYWL,  E,  H,  INIT) 


6-3 .  Arguments 

JE  , 

JM 


PSWL, 

PFWL 


DXWL, 

DYWL 


-  Complex  input  arrays  of  two  elements  each  con¬ 
taining  the  x  and  y  components  of  the  electric 
and  magnetic  surface  current  densities  at  the 
center  of  the  planar  element  as  found  from  E  x  z 
and  z  x  H,  respectively,  where  z  is  the  unit 
normal  to  the  element  and  E,  H  are  the  fields 

at  the  center  of  the  element. 

Real  input  arrays  of  three  elements  each  which 
contain  the  coordinates  Ps(x'/A,  y'/A,  z'/A), 
Pf(x/A,  y/A,  z/A)  of  the  center  of  the  source 
element  and  the  point  at  which  the  field  is  r.o 
be  computed,  respectively. 

-  Real  input  variables  equal  to  the  dimensions  Ax/A, 
Ay/A  of  the  rectangular  source  element. 


I 
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E,  H 


Complex  output  arrays  of  three  elements  each 


containing  the  fields  E,  n  H  at  the  point  P^. 

Note  that  n  H  is  computed  rather  than  H  above 
(to  save  time) . 

INIT  -  Logical  input  variable  which  controls  initiali¬ 

zation  of  various  constants  for  repetitive  calls 
to  the  subroutine:  if  TRUE,  the  constants  are 
computed;  if  FALSE,  the  constants  are  not  com¬ 
puted,  and  their  last  computed  values  are  used. 

6-4 .  Comments  and  Method 

a.  Comment .  The  source  and  field  points  cannot  be  any  closer 
together  than  r  =  .011.  This  restriction  is  necessary  to  prevent  divi¬ 
sion  by  zero  due  to  the  r  ^  variation  of  the  dipole  fields  as  explained 
below.  Actually,  field  points  should  be  removed  to  the  order  of 

/  2  2 

r  =  »' (Ax/X)  +  (Ay/X)  for  validity  of  the  discretized  approximation  to 
the  physical  model. 

b.  Method.  The  subroutine  computation  is  motivated  by  the  problem 

of  computing  the  fields  of  a  rectangular  antenna  aperture  located  in  the 

z=z'  plane  as  illustrated  in  Figure  1.  Let  the  electric  and  magnetic 

fields  E  ,  H  be  specified  at  discrete  points  (x  ,  y  ,  0) .  Then,  at 
— apj  — ap  m  n 

c  .  in 

each  point,  the  equivalent  surface  current  densities  J  and  J  are  given 
by  [1) 


JC  =  z  x  M  -  x  (-H  )  +  y  II  =  x  J*"  +  y  J* 

—  -ap  yap  xap  x  y 


(1) 


m  _  ,  ,  .  in 

J  =  E  x  z  =  x  E  +  y  (-E  )  -  x  J  +  y  J 

-ap  yap  xup  x  y 


(2) 
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The  surface  current  densities  so  defined  can  be  discretized  for 

0 

each  element  AxAy  as  follows.  Consider  the  current  density  .  The 

total  current  entering  the  lower  boundary  and  leaving  the  upper  boundary 
0 

of  the  element  is  J^Ax  and  can  be  regarded  as  an  elementary  dipole  con¬ 
centrated  at  the  center  of  the  element.  The  dipoJe  moment  is 

JS  Ax 

p  =  qi  =  — ^ -  Ay  (3) 

o  30) 


where  q  is  the  charge  and  &  is  the  separation  [2] ,  and  where  the  follow¬ 
ing  relation  for  the  sinusoidal  steady  state  has  been  used: 


I  dt 


I 

ju> 


(4) 


0 

Similar  relations  hold  for  the  other  component  of  J  and,  by  duality, 

for  jm  as  will  be  summarized  below. 

The  next  step  in  the  development  is  to  obtain  expressions  for 
e  m 

the  dipole  fields  of  J  and  J  .  To  facilitate  this  step,  first  con¬ 
sider  the  fields  radiated  by  electric  and  magnetic  dipoles  oriented 

along  the  z  axis  as  shown  in  Table  1  [2].  Note  that  these  expressions 
o 

require  r  in  wavelengths,  and  that  Ax  and  Ay  refer  to  the  element  size 
o 

in  the  original  aperture. 

Matters  are  simplified  if  the  spherical  coordinate  components  of 
Table  1  are  transformed  to  their  corresponding  rectangular  components 
according  to  [3] 
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Table  6—1.  Elementary  Dipole  Fields  of  2— Directed  Currents. 


Electric 


o—  Duality  —  m 


Magnetic 


in  wavelengths 


r  in  wavelengths 


A  =  (A  cos  0  )  sin0  cosiji  +  (A.  sin  0  )  cos  0  cos  <J>  (5a) 

xo  ro  o  o  o  0o  o  o  o 


A  =  (A  cos  0  )  sin8  sin<)>  +  (A-  sin  0  )  cos  0  sin  <p  (5b) 
yo  ro  o  o  o  0o  o  o  o 


A  =  (A  cos  0  )  cos0  -  (A.  sin0  )  sin@  (5c) 

zo  ro  o  o  0o  o  o 


C  =  -  (C,  sin0  )  sin<{>  (5d) 

xo  90  o  o 


C  =  (C.  sin  0  )  cos<{>  (5e) 

yo  <po  o  o 


C  =0  (5f ) 

7.0 

In  the  above,  the  trigonometric  function  in  parentheses  comes  from  the 
field  expressions  in  Table  1;  hence,  the  "minus"  superscript  indicates 
the  field  expression  from  Table  1  without  the  orientation  factor  cos  0 

o 

n  e 

or  sin  0  ,  weighting  e  or  h  ,  and  without  the  current  O'  or  J  . 

o  o  o  zo  zo 

Define  direction  cosines  k  ,  k  ,  k  related  to  0  ,  4>  accord- 

xo  yo  zo  o  o 

ing  to 


k  =  sin  8  cos  $ 
xo  o  o 


k  =  sin  0  sin  41 
yo  o  o 


k  =  cos  0 
zo  o 


Then  Equations  (5)  can  be  rewritten  succinctly  as 


(6a) 


(6b) 


(6c) 
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(7a) 


A 

xo 


(A 


ro 


+  V  k 


xo 


k 

zo 


A  =  (A~  +  a“  )  k  k  (7b) 

yo  ro  6o  yo  zo 


zo 


=  (A  + 
ro 


V 


-  A 
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(7c) 


C  =  -  C  k  (7d) 

xo  <po  yo 


C  =  C,  k  (7e) 

yo  <po  xo 

C  =0  (7f ) 

zo 

Similar  expressions  for  cases  of  x-directed  and  y-directed  di¬ 
poles  may  be  derived  from  those  given  above  merely  redefining  the  axes 
in  Table  1.  When  this  is  done,  the  generalized  expressions  shown  in 
Table  2  result  for  all  three  cases. 

When  both  electric  and  magnetic  currents  are  present  (x-directed 
and  y-directed  components)  the  expressions  for  E  and  H  are  obtained  by 
adding  the  contributions  due  to  each  current  as  given  in  Table  2.  Note 
that  A^,  A(. ,  and  A^_  are  identical  for  both  types  and  directions  of  cur¬ 
rents  so  that  the  expressions  for  the  field  components  may  be  written, 
for  example,  as  follows: 


E 

x 


e 

o 


je((A_ 
x  r 


V 


AJ  +  J  k  k  (A  +  Aft)( 
0  y  x  y  r  0  i 


h  j  tt  n  Jm  k  C 
o  y  z  ij> 


(8a) 
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Table  6-2.  Rectangular  Field  Components  of  Elementary  Dipoles 


Field 

Component 

x-directed 

Dipole  Orientation 
y-directed 

z-directed 

__ 

2 

_ 

_ 

_ 

A 

A  „ 

k 

-  A. 

A  k  k 

A  k  k 

X 

r0 

X 

e 

rt'  x  y 

r0  x  z 

A 

a” 

k 

k 

A  ,  k  -  A. 

A  k  k 

y 

r6 

X 

y 

rf  y  6 

r6  y  z 

A 

A  „ 

k 

k 

A  „  k  k 

-  ,  2 

A  „  k  -A. 

z 

r@ 

X 

z 

r9  y  z 

r0  z  0 

C 

X 

0 

c”  k 

❖  z 

~Ca  k 

4>  y 

C 

y 


k 

z 


0 


k 

x 


C 

z 

where : 


C,  k 


i-0 


<\  *  V 


k 

X 


0 


- 

r  r 
o  o 


-  j2irr 
J,  o 


A 


0 


<2k)2,  -j2’ro 

- )  e 


r 

o 


C 


r 

o 
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Table  6-3.  Fields  of  Elementary  x-I)i  rected  and  y-Direeted  Dipoles 


E  =  Ee  +  Em 


H  =  He  +  Hm 


E  =  e  fJC(A  k2  -  A,)  +  je(k  k  A  .)]  -  h  j'n  k  c 
x  o  x  r.  x  0  y  x  v  re  o  y  z  f 


E  =  e  [Je(k  k  a"  )  +  Je(A  k2  -  A.  )]  +  b  Inn  j”  k  c‘ 
y  o  x  x  y  rO  y  rO  y  h  o  x  z  <J 


E  -  e  [J^Ck  k  A  .)  +  J  (k  k  A  )]  +  h  j-nn  c  ,  (-J  k  +  J  k  ) 

Z  Q  ^  ^  |  <  y  y  2  1  Q  y  ^ 


H  =  e  ^  JG  k  c"  +  h  [  Jm(A  „  k2  -  A  j  +  J1"  k  k  a"  ] 
x  o  n  y  z  >  ol  x  rn  x  e  v  x  y  re 


in  n  —  m  —  rn  —  — 

H  =  -e  J  k  c  +  h  [J  (k  k  A  )  +  J  (A  k"  -  A  )  ] 

y  on  x  z  <p  o  x  x  y  re  y  rn  y 


H  =  e  ^  c  (k  J  -  k  J  )  +  h  |.J  (k  k  A  )  +  J  (k  k  A  )] 
z  o  n  t  y  x  x  y  o  x  x  z  r"  v  y  •»  r- 


,AA, 
(~ )  n 


Where : 


j  (2-0  • 


h  = 
o 


i(2-0' 
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52-53 


54 

55 

56-57 

58-61 

62-65 


67-69 


Compute  E^_  and  due  to  the  x-directed  and  y- 

directed  electric  and  magnetic  currents:  =  JF.  ( 1 )  , 

JS  =  JE (2) ,  Jm  =  JM( 1) ,  Jm'=  JM(2) . 
y  x  y 

2 

Precalculate  (A  +  A_  )  k  -  A„  - 
ro  do  y  do 

See  lines  46-49  above. 


Compute  E  and  H  . 

y  y 

Precalculate  common  variables  for  E  ,  H  . 

z  z 


Compute  E^  and  . 


RETURN 


Error  message  and  halt. 


END 


6-6.  Test  Case 

Selected  test  cases  shown  in  Figure  2-15  of  Reference  1  were  exe¬ 
cuted.  The  square,  4"  x  4",  uniform  aperture  (A=  1.18”)  was  sampled 
at  M=15 ,  N=15  points  in  the  x  and  y  directions,  respectively.  Cases 

were  done  for  E  =  y ( 1 ) ,  H  =0,  and  for  E  =  y(l) ,  H  =  -x(l/n) 

— ap  — ap  — ap  — ap 

In  the  latter  case,  the  amplitudes  obtained  were  twice  as  large  (as 
expected) .  Although  exact  comparison  to  the  graphical  results  in 
Figure  2-15  was  not  possible,  agreement  was  obtained  so  far  as  could  be 
determined.  Some  benchmarks  as  computed  by  Subroutine  DIPOLES  are 
shown  in  Tables  4  through  6. 
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6-8.  Program  Listing.  See  following  panes. 


SUBROUTINE  DIPOLESCJE, JM.PSWL.PFWL.DXWL.DYWL.E.H, INIT)  1 

C  »«»  MODIFIED  1-23-80  TO  INCLUDE  ONLY  1/R  TERMS  ****  NULLIFIED  1-24-80** 

C  SUBR  DIPOLES  COMPUTES  THE  RECTANGULAR  COMPONENTS  OF  THE  FIELDS  E,H  OF  3 

C  ELECTRIC  AND  MAGNETIC  DIPOLES  LOCATED  AT  PSWL(X ' , Y ' , Z ' )  4 

C  AND  ORIENTED  IN  THE  X’  AND  Y'  DIRECTIONS.  THE  FIELDS  ARE  COMPUTED  AT  5 

C  THE  POINT  PFWLCX , Y , Z) .  ALL  DIMENSIONS  ARE  IN  WAVELENGTHS.  MKS  SYSTEM 
C  IS  USED.  FREE  SPACE  (ETA=377  OHMS)  IS  ASSUMED.  7 

COMPLEX  JE(2),JM(2),E(3),H(3),JAY,H0,E0, CPHS , JPI , JAY2PI  8 

COMPLEX  ARO, ATO, CPO, ARTO, CT, ARTOK , CPOK , EOH , HOE  9 

REAL  PSWL ( 3 ) , PFWL ( 3 ) , KX , K Y , KZ  10 

LOGICAL  INIT  11 

C  DXWL, DYWL=X '  AND  Y'  DIMENSIONS  OF  THE  RECTANGULAR  ELEMENT  OVER  WHICH  12 

C  CURRENT  DENSITIES  JE  AND  JM  FLOW  TO  MAKE  THE  DIPOLES.  13 

DATA  TUPI/6. 283185301/. ETA/376. 991 1 185/  14 

IF  (INIT)  GO  TO  1  15 

GO  TO  2  16 

C  COMPUTE  EO.HO  (SEE  DERIVATION  DATED  7-23-79):  17 

1  TUPI2=TUPI**2  18 

JAY= (0 . , 1 . )  19 

E0=DXWL*DYWL*ETA/(JAY«TUPI2)  20 

HO=DXWL*DYWL/( JAY*TUPI2*ETA)  21 

JAY2PI=JAY»TUPI  22 

JPI=JAY*TUPI/2.  23 

HOE=-HO*ETA  24 

EOH=EO/ETA  25 

HO=HO*ETA  26 

E0H=E0H*ETA  27 

C  THE  ABOVE  TWO  LINES  CAUSE  ETA*H  TO  BE  COMPUTED  FOR  USE  IN  RECI.  28 

INIT=. FALSE.  29 

2  X=PFWL( 1 )-PSWL( 1 )  30 

Y =PFWL(2 )-PSWL(2 )  31 

Z=PFWL(3 )-PSWL(3 )  32 

R  =SQRT(X*X+Y*Y+Z*Z)  33 

IF  (R.LT..01)  GO  TO  90  34 

KX=X/R  35 

KY  =Y/R  36 

KZ=Z/R  37 

PHS=AMOD(TUPI*R , TUPI )  38 
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CPHS=CMPLX(0. ,-PHS) 

CPHS=CEXP ( CPHS ) 

ARO=CPHS*( 1 . / R**3+JAY2PI/R**2 ) 

ATO=, 5*( AR0-CPHS*TUPI2/R ) 

CPO=JPI*ARO*R 

ARTO=ARO+ATO 

CT=AR0*KX**2-AT0*( 1 . -KX**2 ) 

ATO=- . 5*CPHS*TUPI2/R 

CPQ=JPI*CPHS*JAY2PI/R 

ARTO=ATO 

CT=-AT0*(1 . -KX*KX) 

ARTOK=ARTO*KX*KY 

CPOK=CPO«KZ 

E(1)=E0*(JE(1) *CT+JE(2 ) *ARTOK )+JM(2 ) *HOE*CPOK 
H( 1 )=H0*( JM( 1 ) *CT+JM(2 ) *ARTOK )+EOH*JE (2) *CPOK 
CT=AR0*KY**2-ATQ*( 1 . -KY**2) 

CT=-AT0*(1.-KY*KY) 

E(2)=E0*(JE( 1 )*ART0K+JE(2)*CT)-JM( 1 )*HOE*CPOK 

H(2)=H0*( JM( 1 )*ART0K+JM(2 )*CT )-JE( 1 )*EOH*CPOK 

ARTOK=ARTO«KY«KZ 

ARTO=ARTO*KX*KZ 

CPOK=CPO*KY 

CPO=CPO»KX 

E(3)=E0*(JE(1 )*ART0+JE(2)*ART0K)+H0E*(JM(  1  )*CPOK 
$~JM(2)*CP0) 

H ( 3 ) =H0* ( JM ( 1 ) * ARTO+JM( 2 ) *ARTOK )+EOH* ( JE ( 1 ) *CPOK 
$-JE(2)*CP0) 

RETURN 

90  WRITE(6,91) 

91  FORMAT ("  ****  R.LT..01  WAVELENGTH  IN  SUBR  DIPOLES— STOP****") 
STOP 

END 


Table  6-4.  Fields  Computed  by  Subroutine  DIPOLES  Along  z-Axis  for  4”  x  4"  Uniform  Aperture 

( A=1 .18") . 
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Table  6-5.  Fields  Computed  by  Subroutine  DIPOLES  Along  x-Axis  at  z=8  inches. 
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Chapter  7 
SUBROUTINE  WALL 

7-1.  Purpose:  To  compute  the  normal  transmission  (T^,T||)  and  reflection 

coefficients  of  a  N-layer  dielectric  sheet  having  thicknesses  d  , 

n 

dielectric  constants  e  ,  and  loss  tangents  tanS  for  each  layer 

rn  n 

when  a  plane  wave  is  incident  at  angle  0. . 

7-2.  Usage:  CALL  WALL  {BETA,  SINE,  D,  ER,  TD,  N,  NN,  TNI,  TN2 , 

RPER,  RPAR) 

7-3.  Arguments 

BETA  -  Real  input  variable  =  2ir/A,  where  A  is  the  free 

space  wavelength. 

SINE  -  Real  input  variable  =  sin  0  . 

D,  -  Real  input  arrays  containing  the  thickness  (cm) , 

ER,  -  dielectric  constant  e^_,  and  loss  tangent  tan6  of 

TD  each  layer. 

N  -  Integer  input  variable  equal  to  the  number  c  f 

layers . 

NN  -  Integer  input  =  N+l. 

TNI ,TN2  -  Complex  output  variables  equal  to  the  normal 

voltage  transmission  coefficients  for  the  components 
of  the  incident  electric  field  perpendicular  to  and 
parallel  to  the  plane  of  incidence,  respectively. 

RPER, RPAR  -  Complex  output  variables  equal  to  the  reflection 

coefficients  R  ,  R,.  . 

i  II 
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7-4 .  Comment  and  Method 


a.  Layer  1  is  the  first  layer  on  the  exit  side  of  the  panel;  layer 
N  is  the  first  layer  on  the  incident  side.  T^ ,  T  j  |  have  the  same  value 
for  either  side  of  the  panel  being  the  incident  side;  however,  R^,  R||  are 
different  (in  phase)  for  the  two  cases. 

b.  The  details  of  the  method  are  presented  in  Appendix  E  of  Refer¬ 
ence  1 . 


7-5.  Program  Flow:  See  Reference  1. 

7-6.  Test  Case:  None. 

7-7.  References 

1.  G.  K.  Huddleston,  "Radome  Analysis  Computer  Program:  Ray  Tracing 
Formulation",  Technical  Report  for  JHU/APL,  Contract  No.  601053, 
November  1979. 

7-8.  Program  Listing:  See  following  pages. 


on  oooooooooooo 


X 


SUBROUTINE  WALL (BETA , SINE , D, ER , TD, N, NN , TN 1 ,TN2, RPER , RPAR )  1 

C  SUBROUTINE  WALL  COMPUTES  THE  TRANSMISSION  AND  REFLECTION  2 

C  COEFFICIENTS  FOR  AN  N  LAYER,  PLANE  DIELECTRIC  PANEL  FOR  PLANE  3 

C  WAVE  INCIDENT  AT  SINE(ANGLE)  FOR  PERPENDICULAR  AND  4 

C  PARALLEL  POLARIZATIONS.  5 

C  PARAMETERS  OF  THE  WALL:  N=  THE  NUMBER  OF  LAYERS  6 

C  NN=  N+1  REQUIRED  TO  DIMENSION  ARRAYS  7 

C  D=  THICKNESS  OF  EACH  LAYER  IN  CENTIMETERS  8 

ER  =  RELATIVE  DIELECTRIC  CONSTANT  OF  EACH  LAY  9 

TD=  THE  LOSS  TANGENT  FOR  EACH  LAYER  10 

TN1.TN2  ARE  THE  NORMAL  VOLTAGE  XMN  COEFFICIENTS;  TPER.TPAR  ARE  THE  1 1 

INSERTION  VOLTAGE  TRANSMISSION  COEFFICIENTS.  IT  IS  IMPORTANT  TO  12 

NOTE  THAT  THE  XMN  COEFS  ARE  THE  SAME  FOR  PLANE  WAVE  INCIDENT  FROM  13 

EITHER  SIDE  OF  THE  STRATIFIED  DIELECTRIC  PANEL  IMMERSED  IN  FREE  SPACE;  14 

HOWEVER,  THE  REFLECTION  COEFS  ARE  NOT.  THAT  IS,  FOR  COMPUTING  RPER,  15 

RPAR,  THE  ORDERING  OF  ER(NN) ,TD(NN)  IS  IMPORTANT  WITH  LAYER  1  BEING  16 

THE  FIRST  LAYER  ON  THE  EXIT  SIDE,  LAYER  N  BEING  THE  FIRST  LAYER  ON  THE  17 

INCIDENT  SIDE.  LAYER  NN  AND  LAYER  0  ARE  JUST  FREE  SPACE  LAYERS  18 

OF  SEMI-INFINITE  DEPTH.  19 

E,G,R1,R2,  ARE  ARRAYS  USED  IN  THE  SUBROUTINE  HAVING  NN  DIM"L  LIMITS  20 

COMPLEX  E(6),G(6),R1(6),R2(6),GG,EE,RR1,RR2,AA1, AA2, XI , X2,  21 

$X3.X4,Y1,Y2,Y3,Y4,U1,U2,U3,U4,V1,V2,V3,V4,P1,P2,P3,P4,Q1,Q2,Q3.Q4  22 

COMPLEX  TPER.TPAR, RPER, RPAR, U,V, TNI, TN2  23 

DIMENSION  ER(NN),TD(NN),D(N)  24 

ER(NN)= 1 .0  25 

TD(NN)=0.  26 

DO  50  1=1, NN  27 

50  E(I)=CMPLX(ER(I),-ER(I)«TD(I))  28 

AB=BETA*0. 70707070707071  29 

30 

CALCULATE  TOTAL  THICKNESS  OF  WALL  IN  CM  31 

C  32 

DTOTAL=0 . 0  33 

DO  200  1=1, N  34 

200  DTOTAL=DTOTAL+D ( I )  35 

C  S  IS  THE  SINE  OF  THE  ANGLE  SQUARED  36 

C  C  IS  THE  COSINE  OF  THE  ANGLE  37 

S=SINE*SINE  38 
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C=SQRT ( 1 .0-S) 

AD=ER(1 ,-S 
ET=ER(1 )*TD(1 ) 

SR=SQRT(AD*AD+ET*ET) 

IF(SR-AD)  76,76,77 

76  A=0. 

GO  TO  78 

77  A=AB*SQRT(SR-AD) 

78  B=AB*SQRT(SR+AD) 

G ( 1 ) =CMPLX( A, B) 

GG=CMPLX(0 . 0, BETA*C) 

EE= 1 . 0 

SUM=0. 

SUM=SUM+D(1 )/SQRT(AD) 

RR 1  =  (G ( 1 )-GG)/(G(1 )+GG) 

RR2= (EE*G ( 1 )-E ( 1 ) *GG) /(EE*G( 1 )+E ( 1 ) *GG) 

DO  84  1=1, N 
11=1+1 
AD=ER(II)-S 
ET=ER(II)*TD(II) 

IF  (I-N)  176,177,177 

176  SuM=3UM+D (II ) /SQRT ( AD ) 

177  CONTINUE 
SR=SQRT(AD*AD+ET*ET) 

IF(SR-AD)  79.79.80 

79  A=0. 

GO  TO  81 

80  A=AB*SQRT(SR-AD) 

81  B=AB*SQRT(SR+AD) 

G ( I I ) =CMPLX ( A , B) 

R1(I)  =  (G(II)-G(I))/(G(II)+G(D) 

84  R2(I)=(E(I)*G(II)-ECII)*G(I))/(E(I)*G(II)+E(II)*G(I)) 
SUM=S«SUM 

AA 1=1. O-RR 1 
AA2=1 . 0-RR2 
DO  85  1=1. N 
AA1=AA1*(1 ,0-RKI)) 

85  AA2=AA2*(1 ,0-R2(I) ) 


39 

no 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 
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AA 1  =  1 . 0/  AA1 

77 

AA2=1 . 0/  AA2 

78 

U=-G( 1 )*D( 1 ) 

79 

V  =G  ( 1 )  *D  ( 1 ) 

80 

XI =CEXP(U ) 

81 

X4=CEXP (V ) 

82 

X2=-RR1»X4 

83 

X3=-RR1*X1 

84 

Y 1  =X  1 

85 

Y4=X4 

86 

Y2=-RR2*Y4 

87 

Y3=-RR2*Y 1 

88 

DO  105  1=2 , NN 

89 

IF(I-NN)  95,90,1 

90 

Ulsl.O 

91 

U2=-R1 (N ) 

92 

U3=-R1(N) 

93 

U4=1 .0 

94 

VIsl.O 

95 

V2=-R2(N) 

96 

V3=-R2(N) 

97 

V4=1 . 0 

98 

GO  TO  100 

99 

11=1-1 

100 

U=-G(I )*D (I ) 

101 

V=G(I)*D(I) 

102 

U1=CEXP(U) 

103 

U4=CEXP (V ) 

104 

U2=-R1 (II)*U4 

105 

U3--R1 (II)*U1 

106 

VI  =U  1 

107 

V4=U4 

108 

V2=-R2(II)*V4 

109 

V3=-R2(II )*V 1 

110 

P1=X1»U1+X2*U3 

111 

P2=X1*U2+X2*U4 

112 

P3=X3*U1+X4»U3 

113 

P4=X3*U2+X4»U4 

114 
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Q 1 =Y 1 *V  1  +Y2*V  3  115 

Q2=Y1»V2+Y2*V4  116 

Q3=Y3*V1+Y4«V3  117 

Q4=Y3*V2+Y4*V4  118 

X1=P1  119 

X2=P2  120 

X3=P3  121 

X4=P4  122 

Y1=Q1  123 

Y2=Q2  124 

Y5=Q3  125 

105  Y4=Q4  126 

RPER=-X3/X4  127 

TN1.TN2  ARE  NORMAL  VOLTAGE  XMN  COEFFICIENTS.  128 

RPAR=-Y3/Y4  129 

TN 1  =  (X 1+X2*RPER  ) *AA  1  130 

U :CMPLX (0 . 0 , -SUM* BETA )  131 

U=CEXP (U )  132 

TPER.TPAR  HERE  ARE  VOLTAGE  XMN  COEFFICIENTS  AT  EXIT  POINT  OF  RAY.  133 

TPER  =TN  1  *1)  134 

TN2=(Y 1+Y2*RPAR)*AA2  135 

TPAR=TN2*U  136 

MODIFY  TRANSMISSION  COEFFICIENTS  FOR  INSERTION  137 

U=CMPLX(O.O.BETA»DTOTAL*C)  138 

U=CEXP(U)  139 

TPER=TN1*U  140 

TPAR=TN2*U  141 

1  CONTINUE  142 

300  RETURN  143 

END  144 


Chapter  8 
SUBROUTINE  CAXCB 


8-1.  Purpose:  To  compute  the  complex  vector  cross  product  C  =  A  x  B, 

where  A  and  B  are  a  complex  vectors  expressed  in  rectangular  coor¬ 
dinates  . 

8-2.  Usage:  CALL  CAXCB  (A,  B,  C) 

8-3.  Arguments 


Complex  input  array  containing  the  rectangular 
components  of  the  vector  A  =  x  A  +yA  +  z  A  ; 

*  y  = 

i .e. ,  A  (A  ,  A  ,  A  ) . 

x  y  z 

Complex  input  array  B  (B^,  B^,  B^)  representing 
the  vector  B. 

Complex  output  array  C  (C^ ,  C  ,  C^)  representing 


the  vector  C  =  A  x  li. 
8-4.  Comment  and  Method:  None 
8-3.  Program  Flow:  See  listing  below. 

8-6.  Test  Case:  None. 

8-7.  References:  None. 

8-8.  Program  Listing:  See  following  page. 
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SUBROUTINE  CAXCB(A,B,C)  1 

COMPLEX  A(3),B(3),C(3)  2 

SUBR  CAXCB  COMPUTES  THE  VECTOR  CROSS  PRODUCT  C=AXB  OF  3 

TWO  COMPLEX  VECTORS  A  AND  B  4 

C(1 )=A(2)»B(3)-A(3)*B(2)  5 

C(2)=A(3)*B(1)-A(1)»B(3)  6 

C ( 3 ) =A ( 1 )*B(2)-A(2)*BO )  7 

RETURN  8 

END  9 


